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Abstract 

The number of nowhere-zero Zq flows on a graph G can be shown to be a polynomial 
in Q, defining the flow polynomial $g{Q)- According to Tutte's five-flow conjecture, 
$g(5) > for any bridgeless G. A conjecture by Welsh that <&g(Q) has no real roots 
for Q 6 (4, oo) was recently disproved by Haggard, Pearce and Royle. These authors 
conjectured the absence of roots for Q £ [5, oo). We study the real and complex 
roots of <&g(Q) f° r a family of non-planar cubic graphs known as generalised Petersen 
graphs G(n, k). We show that the modified conjecture on real flow roots is also false, by 
exhibiting infinitely many real flow roots Q > 5 within the class G(nk, k). In particular, 
we compute explicitly G(119, 7) showing that it has real roots at Q ~ 5.0000197675 and 
Q ~ 5.1653424423. We moreover prove that the graph families G(6n, 6) and G(7n, 7) 
possess real flow roots that accumulate at Q = 5 as n — > oo (in the latter case from 
above and below); and that Q c (7) ~ 5.2352605291 is an accumulation point of real 
zeros of the flow polynomials for G(7n,7) as n — > oo. 



Key Words: Nowhere-zero flows; flow polynomial; flow roots; Tutte's five-flow conjecture; 
Petersen graph; transfer matrix. 



1 Introduction 



The four-colour problem is probably the best known problem of graph theory, and has 
fascinated generations of mathematicians and laymen since it was first stated in a letter 
from Augustus de Morgan to Sir William Rowan Hamilton dated 23 October 1852. The 
suspense was resolved in late July 1976 when Appel and Haken [1] announced the proof of 
the four-colour theorem: Every planar graph admits a 4-vertex-colouring. 

A quantitative version of the problem was proposed by Birkhoff [14] in 1912: Given an 
arbitrary graph G and a set of Q colours, how many proper vertex colourings does G admit? 
This defines the chromatic polynomial Xg(Q)- The four-colour theorem can then be stated: 
If G is planar, xg(4) > 0. 

Although initially defined for Q G N, it is easily proved [14,58,59] that Xg(Q) is i n f &c t a 
polynomial in Q, and can as such be evaluated for any Q G C. This suggests an algebraic or 
even analytic approach to the colouring problem. There exist many studies of the location 
in C of the roots of Xg(Q), henceforth called chromatic roots. These studies concern either 
specific graphs, or all planar graphs, or other infinite families of graphs. 

Birkhoff and Lewis [15] have made the following conjecture: If G is planar, Xg(Q) > 
for Q G [4, oo). Obviously, the statement of this conjecture is stronger than the four-colour 
theorem, but unfortunately has not yet been turned into a theorem. The corresponding 
result for Q G [5, oo) has however been proved by the same authors [15] (see also [37,50,60]). 

Beraha and Kahane [8] have exhibited an infinite family of planar graphs for which Q = 4 
can be proved to be an accumulation point of complex chromatic roots. In that sense, the 
four-colour theorem is "almost false". Improving on this, Royle [41] has proved, for a slightly 
different family, that Q = 4 is also an accumulation point of real chromatic roots (converging 
to Q = 4 from below). Finally, Sokal [48] has proved that for a specific family of planar 
graphs (generalised ©-graphs) chromatic roots are dense in C (except perhaps in the disc 
|Q-1|<1). 

A close cousin of the chromatic polynomial is the so-called flow polynomial <&g(Q)- Let 
G = (V, E) be an arbitrary (not necessarily planar) graph G with vertex set V and edge 
set E, and let T be an additive Abelian group. A T-flow on G is a map 0: E \-> T that 
attributes a variable 0(e) to each edge e G E, subject to the conservation of these variables 
at each vertex, with respect to an arbitrary chosen orientation of E. An elementary example 
of a flow is the current in an electrical network, in which case the conservation constraint 
is known as Kirchhoff's first law [34]. 

A nowhere zero T-flow is a T-flow such that 0(e) ^ for all e G E [33,54,63]. If 
T is a finite Abelian group of order Q, it can be shown that the number of nowhere zero 
T-flows depends only on Q (not on the specific structure of the group T), and it is in fact 
the restriction to Q G N of a polynomial in Q called the flow polynomial <&g{Q) [52]. One 
can then again extend the definition to Q G C and study the location of (real or complex) 
flow roots. 

A nowhere zero Q-flow of G is a nowhere zero Z-Aow such that |0(e)| < \Q — 1| for all 
e G E. Tutte [52] showed that G has a nowhere zero Q-flow if and only if it has a nowhere 
zero Zq-Aow; but these two concepts are different! Tutte's result immediately implies the 
following interesting (but far from obvious) property of nowhere- zero Zq-Aows [53] : 

Proposition 1.1 If$ G (Q) > for some Q G N, then $g(Q') > for all integers Q' > Q. 
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When G is planar, one has [51] the duality relation \g* (Q) = Q $g(Q), where G* denotes 
the dual graph. The properties of &g(Q) thus follow from those of Xg*(Q)- 

It is worth stressing that the Birkhoff-Lewis theorem [15] provides a uniform upper bound 
for the real zeros of the chromatic polynomial of all planar graphs, namely Q = 5. However, 
such an upper bound (if it actually exists at all!) is not known for the real zeros of the flow 
polynomial of arbitrary bridgeless graphs. 1 The existence of such uniform upper bound and 
its value, if it does exist, are long-standing open problems in Combinatorics. 

Consider now arbitrary (not necessarily planar) bridgeless graphs G. Because there exist 
graphs not admitting a nowhere zero 4-flow, 2 the strongest possible results for integer and 
real flow roots are given respectively by the following two well-known conjectures: 

Conjecture 1.2 (Tutte's five-flow conjecture [33,53,54]) For any bridgeless graph G, 
$g(5)>0. 

Remark. This conjecture implies that <&g(Q) > for all integers Q > 5 by Proposition 1.1. 
Conjecture 1.3 (Welsh [57]) For any bridgeless graph G, &g{Q) > for Q e (4, oo). 

It should be noted that the Welsh conjecture parallels that of Birkhoff and Lewis for 
the chromatic polynomial. It does however not include the case Q — 4. Some results by 
Jackson on zero-free intervals for the flow polynomials of cubic graphs [24, 25] also suggest 
this close parallelism between xg(Q) f° r planar G and $g(<3) for arbitrary G. (Note that 
both polynomials are evaluated at the same value of Q.) 

A number of weaker results have been proved over the years, notably: 

Theorem 1.4 (Seymour [46]) For any bridgeless graph G, &g(6) > 0. 

Theorem 1.5 (Steinberg [49]) For any bridgeless graph G that is embeddable in the pro- 
jective plane, $g(5) > 0. 

An immediate Corollary of Seymour's theorem (using Proposition 1.1) is that <&g{Q) > 
for all integers Q > 6. Thus Q = 5 is a uniform upper bound for integer flow roots. But the 
above results give no clue about the existence of a uniform upper bound for real flow roots. 

A first step into proving (or disproving) Conjectures 1.2-1.3 consists in studying the flow 
roots of "small" graphs. By computing the flow roots of small graphs with high girth (up 
to 32 vertices and girth at least 7), 3 Haggard, Pearce, and Royle [22] have very recently 
found an explicit counterexample to the Welsh conjecture: the flow polynomial of the gen- 
eralised Petersen graph G(16,6) has two real roots larger than Q = 4: <5i ~ 4.0252205, 
and Q 2 ~ 4.2331455. However, the same authors conjectured the following modification of 
Conjecture 1.3, in which 4 is replaced by 5, and the endpoint Q = 5 is now included in 
accordance with Tutte's five-flow conjecture: 

1 Obviously, $g(<5) = if G has a bridge, because of the "nowhere zero" condition. 

2 Indeed, the Petersen graph — which is a special case G(5, 2) of the generalised Petersen graphs G(n, k) to 
be defined in section 3 below — has the flow polynomial $g(5.2)(Q) = (<5-l)(Q-2)(Q-3)(Q-4)(Q 2 -5(5+10), 
which vanishes at Q = 4. 

3 Note that Kochol [35] proved that the smallest counterexample of Tutte's five-flow conjecture should 
have girth at least 9. 
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Conjecture 1.6 ( Haggard-Pear ce-Royle [22]) For any bridgeless graph G, $g(<3) > 
for Q E [5, oo). 

Even though the naive parallelism between Xg{Q) an d &g(Q) has been invalidated by 
the above-mentioned counterexample to Conjecture 1.3, a related line of reasoning would be 
that all these conjectures and theorems might be related by replacing Q for the chromatic 
polynomial with Q + l for the flow polynomial. Thus, the four-colour theorem [1] "translates" 
into the Tutte five-flow conjecture [53,54]. Similarly, Conjecture 1.6 is the "translated" 
version of the Birkhoff-Lewis conjecture. (Note that the translation of Royle's result [41] is 
consistent with the fact that Conjecture 1.3 [57] is false.) 

In this paper we study the flow polynomial on the infinite family of graphs known as the 
generalised Petersen graphs. Our main results are the following: 

Theorem 1.7 The value Q = 5 is an isolated accumulation point of real zeros of the flow 
polynomial $g(Q) for the families of bridgeless graphs G(6n,6) and G(7n,7) with n > 3. 
Moreover, 

(a) There is a sequence of real zeros {Q n } of the flow polynomials $G(6n,6) that converges 
to Q = 5 from below. 

(b) There is a sequence of real zeros {Q n } of the flow polynomials $G(7n,7) that converges 
to Q = 5. The sub-sequence with odd (resp. even) n converges to Q = 5 from above 
(resp. below). 

Theorem 1.8 

(a) The bridgeless graph G(119, 7) has flow roots at Q 5.00002 and Q 5.16534 (where 
ps means "within 10~ 5 "). 

(b) The value Q c (7) ~ 5.235261 (where ~ means "within 10~ 6 ") is an accumulation point 
of real zeros of the flow polynomials $G(7n,7) (Q) ■ ^ n particular, the sub-sequence for 
odd n of the real zeros {Q n } of the flow polynomials $G(7n,7) converges to Q C (J) from 
below. 

Remark. The largest real flow root we have explicitly found is Q pa 5.1653424423. 

Thus, the Welsh conjecture and the "translated Birkhoff-Lewis conjecture" are false, and 
the Tutte conjecture is "almost false" (in the same sense that the four-colour theorem is 
"almost false" [8,41]). On the other hand, Theorem 1.7 includes the "translated version" of 
that of Royle [41] for real chromatic roots. 

In this work, we have considered the family of graphs G(nk, k) with k < 7 and n > 2. For 
each k, we have located the set of accumulation points in the complex Q-plane of the roots 
of the flow polynomial 3>G(nfc,fc)> as n ~ * 00 • Most accumulation points belong to limiting 
curves B^; and in particular, we are interested in locating the points Q c (k), defined as the 
largest real value where the limiting curves Bk cross the real axis. (These points are likely to 
be accumulation points as n — > oo of real zeros, as in Theorem 1.8(b); but not always: see 
section 6.3.) We have been able to obtain the values of Q c (k) for k < 11; and the numerical 
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extrapolation of these values to k — > oo yields Hindoo Q c (k) = Qo ~ 5.75. We expect that 
this value is the largest real accumulation point one can get from the family G(nk, k). 

Based on this — and on the failure of Conjectures 1.3 and 1.6 — we venture the following 
weaker conjecture: 

Conjecture 1.9 For any bridgeless graph G, &g(Q) > for Q e [6, oo). 

This translated version of the Birkhoff-Lewis theorem [15] bears the same relation to 
the Seymour theorem [46] as the original Birkhoff-Lewis theorem bears to the 5-colour 
theorem, or the Birkhoff-Lewis conjecture does to the 4-colour theorem; this can be taken 
as (weak) evidence that Conjecture 1.9 might be true. On the other hand, firm believers 
of the translation mechanism outlined above would conclude — by reverse translation of our 
Theorem 1.8 (and its extension) — that the Birkhoff-Lewis conjecture is false! We do not, 
however, possess any evidence whatsoever in support of this conclusion, based on actual 
studies of the chromatic polynomial. 

Note that we do in fact exhibit infinitely many flow roots in an interval Q G [4, Q ] with 
Qo > 5. This means, loosely speaking, that one should look for a purely combinatorial 
proof of Tutte's 5-flow conjecture. This is exactly the same situation as for the 4- and 5- 
colour theorems; they hold true even though one can find families of graphs with real roots 
approaching Q = 4 from below [41], and other families with complex roots approaching 
densely to Q = 4 and Q = 5 [48]. If Conjecture 1.9 turns out to be false, then it is very 
plausible that there exist no upper bound for real flow roots of arbitrary bridgeless graphs. 

On a more technical level, we exhibit a method for computing exactly the flow polynomial 
on very large generalised Petersen graphs (which can readily be adapted to other similar 
graph families). This method relies on a transfer matrix construction similar to the one 
employed in our previous work [26,27] on the chromatic polynomial for graphs with periodic 
longitudinal boundary conditions. 

The paper is organised as follows. In section 2 we define the flow polynomial carefully 
and exhibit its relation to the Q-state Potts model. Building on this, we show in section 3 
how the flow polynomial on generalised Petersen graphs can be built by a transfer matrix 
construction. Our results, given in section 4, are obtained by implementing this construction 
on a computer and pushing the computation to as large graphs as possible. Note that 
although obtained by computational means, the flow polynomials are exact and involve no 
approximation whatsoever. In section 5 we introduce the Beraha-Kahane-Weiss theorem, 
that plays an important role in establishing our results. In section 6, we describe our analytic 
findings about the real zeros of the flow polynomial for this family of graphs, and in section 7, 
we summarise the more experimental (i.e., numerical) aspects of our work. To conclude, 
section 8 contains a discussion of our results. In Appendix A, we prove some technical 
lemmas included in the text, and in Appendix B, we give the coefficients of the generalised 
Petersen graph P(119, 7). 

2 Flow polynomial 

Let G = (V, E) be a connected graph and T be an Abelian group. Assign an arbitrary 
orientation to each edge e e E. With respect to any fixed vertex i e V, the edges Ei incident 
on i can then be characterised as either ingoing or outgoing: E^ = E™ U E° ut . 
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A T-flow on G is a map 0: E h-» T that attributes a variable 0(e) to each edge e E E, 
subject to the constraint 

E <^ e ) = E <^ e ) (2- 1 ) 

e£E\ n e£E° ut 

for any i e V. The edge orientation is actually immaterial in these definitions: if one wants 
to change the orientation of an edge eo, it suffices to change simultaneously the sign of the 
flow along that edge, 0(e o ) — > — 0(e o ). 

A nowhere zero T-flow is a T-flow such that 0(e) ^ for all e e E. If T is a /mite 
Abelian group of order Q, then the number of nowhere zero T-flows depends only on Q, and 
is called the flow polynomial <&g(Q), for reasons we shall exhibit in the next paragraph. In 
particular, a Zq-flow (resp. a nowhere zero Zq-flow) on G is a map : E 1 >-)■ {0, 1, . . . , Q — 1} 
(resp. 0: t-> {1, 2, . . . , Q — 1}) for which the constraint (2.1) is imposed modulo Q. 

Clearly, the total number of Zg-flows on G is Q°^ E \ where for any subset E' C E 1 , 
c(£") denotes the number of independent cycles (cyclomatic number) in the induced graph 
G' = (V,E'). To obtain the number of nowhere zero Zq-Aows, we first subtract for each 
e G E the flows for which 0(e) = 0. Since flows with two zero-flow edges will be subtracted 
off twice, these must be put back in the sum, and proceeding by inclusion-exclusion we 
find [61] 

$ G (Q) = J2(-l) lEHE ' l Q c(E,) - (2.2) 

E'CE 

By this result, <&g(Q) 1S indeed a polynomial in Q, and (2.2) is used to extend the definition 
to Q E C. 

Meanwhile, recall the partition function of the Q-state Potts model [39] 

Z G (Q,v) = E II e K <^ U) ) , (2.3) 
o- (ij)eE 

where the map a: V — > {0, 1, . . . , Q — 1} is called the spin, and K is the coupling constant. 
The Kronecker delta function 5(x,y) is defined by S(x,y) = 1 if x = y, and S(x,y) = 
otherwise. We have introduced the convenient parameter v = e K — 1. By expanding the 
edge product and performing the sum over a, one recovers the partition function in the 
Fortuin-Kasteleyn cluster representation [20] 

Zg(Q,v) = E^V (E ' } , (2.4) 

E'CE 

where k(E') is the number of connected components in G' = (V, E'). 

Graph theorists will recognise in (2.4) [a reparametrisation of] the Tutte polynomial [53] 
and interpret (rasa vertex colouring. Proper vertex colourings, i.e., those for which adjacent 
vertices are coloured differently, are obtained for K — > — oo, and therefore 

Xg(Q) = Z G (Q,-1) (2.5) 

is the chromatic polynomial. 

Setting instead v = —Q in (2.4), and using the topological identity 

k(E') = \V\ - \E'\ +c(E'), (2.6) 
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one establishes the connection with the flow polynomial 

<MQ) = (-l) lEl Q- lvl Z G (Q,-Q). (2.7) 

In the case where G is planar, let G* denote the dual graph. Recall the fundamental 
duality relation [62] of the Potts model partition function 

Z G (Q,v) = KZ G *(Q,v*), (2.8) 

where v* is the dual of v 

vv* = Q, (2.9) 

and the proportionality factor is 

K = QHVI^I = q\v\-\e\-i v \e\ (21Q) 

Noticing that v = —Q is dual to v — — 1 by (2.9) furnishes a relation between the flow 
polynomial of G and the chromatic polynomial of G*. Indeed, using (2.8) and (2.10) we 
have [51] 

xg*(Q) = Q$g{Q). (2.H) 

Note that <&g(Q) = if G contains a bridge e € E. Indeed, by the constraint (2.1) one 
would have 0(e o ) = 0, preventing the existence of a nowhere zero flow. This observation is 
consistent with the duality relation (2.11), since a bridge in G implies a loop in G*. 

Alternatively, the relation (2.11) can be proved by noting that there exists an obvious 
bijection between the nowhere zero Zg-flows on G and the proper colourings of the faces of 
G, with the colour on one face being fixed. Indeed, let be a flow on G. Then, turning 
around a vertex, each time one moves from a face i to an adjacent face j, if the separating 
edge e is seen oriented to the right (resp. left), its flow variable 0(e) defines the colour 
difference Cj — Ci = 0(e) (resp. Cj — Cj = —0(e)). Starting from the face with fixed colour, 
these differences define the face colouring of the whole graph. The mapping from proper 
colourings to flows follows similarly. 

It is useful to note that for any bridgeless 3-connected graph G, one can deduce from 
(2.4)/(2.7) that <&g(Q) is a polynomial in Q of degree \E\ — \V\ + 1 in Q, and that the first 
two coefficients of &g(Q) ar e given by 

MQ) = Q lEHVl+l - \E\ Q' £ '-I y l + . . . (2.12) 

The first term comes from the fact that there is a unique spanning graph (V, E') in (2.4) / (2.7) 
with E' = E. The second term is given by the contribution of the \E\ spanning subgraphs 
(V, E') with E' = E \ e for each e e E, and the observation that k(E') = 1 since G 
is connected and bridgeless. Notice that if we consider the spanning graph (V, E') with 
E 1 \ {e, e'} for any two distinct edges e, e' G E, the 3-connectedness of G guarantees that 
the next term in (2.12) with be of order Q'^' - ^' -1 , as in this case we also have k(E') = 1. 
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3 Transfer matrix for flow polynomials of generalised 
Petersen graphs 



3.1 Generalised Petersen graphs 

The goal of this paper is to evaluate the flow polynomial on a family of graphs G(n, k) 
called generalised Petersen graphs and defined as follows: Let n, k be positive integers such 
that n > k. Then G(n, k) is a cubic graph with 2n vertices denoted i p and j p for p = 
1,2,..., n: i.e., 

V(G(n, k)) = {i 1 ,...,i n ,j 1 ,...,j n }. (3.1) 

The edge set consists of 3n edges (i p j p ), (i p i p+ i), (j p j p +k), f° r P — 1,2, . . ■ ,n, and with all 
indices considered modulo n: i.e., 

E(G(n, k)) = {(ip,j P ), (ipip+i), (j P j P +k) I 1 < P < n} . (3.2) 

Note that G(n, k) is simple for n ^ 2k; but it has double edges when n = 2k. These graphs 
were introduced by Watkins [56]. As an example, (2(12,4) can be drawn as follows: 



The graphs G(n, k) are clearly bridgeless. They are non-planar for all pairs (n, k), except for 
the case (3, 2), and the two sub-families (p, 1) and (2p, 2) with p > 1. They have girth 8 for 
n and k sufficiently large. We have thus a two-parameter family of non-planar cubic graphs 
with high girth, and based on exhaustive studies of small graphs [22] we expect this family 
to produce large real flow roots. 

Remark. Note that we do not expect that any of these roots to be Q = 5. In fact, the 
generalised Petersen graphs satisfy the following 

Lemma 3.1 For every generalised Petersen graph G{n,k) with n,k positive integers such 
that k < n, $G(n,fe)(5) > 0. In fact, every graph G(n, k) other than the ordinary G(5, 2) has 

$G(n,fc)(4) > 0. 

Proof. It is well known [16] that every generalised Petersen graph G(n, k) (with the excep- 
tion of the Petersen graph (7(5,2) itself), admits a Tait colouring: i.e., an edge 3-colouring 
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such that at every vertex, the three incident edges take distinct colours. 4 The existence of 
such edge 3-colourings is equivalent, for cubic and loop less graphs, to the existence of a 
nowhere-zero 4-flow for the same graph [32, Proposition 2(b)]. Therefore, for all G(n,k), 
except the Petersen graph G(5, 2), <J> G ( n fc )(4) > 0, and furthermore, $G( n ,fc)(5) > by Propo- 
sition 1.1. The case (2(5,2) is dealt with directly: from the exact expression for $g(5,2) (see 
footnote 2), we conclude that $g(5,2)(5) = 240 > 0. I 

We shall however show that the five-flow conjecture is "almost false", in the sense of 
Theorem 1.7. 

3.2 Potts model transfer matrix 

We wish to evaluate Z G ^ n ^(Q,v) — of which the flow polynomial &G(n,k)(Q) is a special 
case — by a transfer matrix construction. 

Contrary to an often repeated but false statement, evaluating Z G (Q,v) by a transfer 
matrix construction is possible for any graph G, and does not require G to consist of a 
number of identical layers [7]. However, when G does have a layered structure — as is the 
case here — Zq(Q,v) can be computed by the repeated application of the same transfer 
matrix. 

Supposing for simplicity that n is a multiple of k, the generalised Petersen graph G(n, k) 
can be redrawn as in Figure 1. This turns G{n,k) into a graph of n/k identical layers of 
width L = k + 1 vertices with periodic boundary conditions in the vertical direction. (We 
shall henceforth refer to this as periodic longitudinal boundary conditions, in accordance 
with the fact that the transfer matrix builds up the graph vertically.) We now claim that 
this implies that Za( n ,k)(Q,v) can be written as a Markov trace 

Z G (n,k)(Q,v) = Tr(T L ) n/k (3.3) 

of the (n/k)th power of a transfer matrix Tl to be defined shortly. 

In general, for a layered graph of width L, T L acts on basis states A L which are set 
partitions of 2L points {1, 2, . . . , L, 1', 2', . . . , L'}. These basis states can be depicted as 
partition diagrams, which are hypergraphs on 2L vertices, drawn inside a rectangle with L 
vertices (labelled 1', 2', . . . , V) on top and L vertices (labelled 1, 2, . . . , L) on bottom. Each 
hyperedge represents one block in the partition. A block that contains at least one vertex 
from both the top and bottom rows is called a link. A block containing precisely one vertex 
is called a singleton. The number of links in a diagram d is denoted £(d). 



4 It is worth noting that the definition of the generalised Petersen graph G(n, k) in Refs. [16,56] explicitly 
excludes the case n = 2k. However, it is easy to see that any G(2k, k) has a Tait colouring: e.g., the edges 
(ip>jp) take colour 1, the edges (i p , i p +i) take alternatively colours 2 and 3 (as p goes from 1 ton), and for 
each p, one of the double edges (j p ,j p +k) takes colour 2, and the other edge, colour 3. 
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(a) (b) 

Figure 1: (a) Generalised Petersen graph G(n, k), here with k — 4. There are n layers of two 
vertices in the vertical direction, but some edges link layers at distance k. The boundary 
conditions are periodic in the vertical direction, (b) When n G kN, G(n, k) can be redrawn 
as shown. There are nj k layers of width k + 1 vertices, each comprising a total of 2k vertices. 
All edges now link vertices within the same layer, or in two adjacent layers. 



The following example with L = 4 



1' 2' 3' 4 



d\ = 




represents the partition (12'3')(231'4')(4). It has two links and one singleton. 

The multiplication d = d 2 -d\ of two partition diagrams is defined by stacking the diagrams 
vertically. Specifically, the top row of d 2 becomes the top row of d, the bottom row of d\ 
becomes the bottom row of d, and the top row of d\ is identified with the bottom row of d 2 . 
Any blocks not containing points in the top or bottom rows of d are removed in the process. 
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This gives, for example: 

V 2' 3' 4' 1" 2" 3" 4" V 2' 3' 4' 




This diagram multiplication turns Al into an associative partition monoid [23, 36] with 
identity I = (11') (22') • • • (LL'). Observe that 

£(d 2 -d 1 ) < min(e(d 2 ),e(d 1 )) . (3.4) 

The idea is now that these diagrams will represent the edge subset appearing in the 
cluster representation (2.4) of the Potts model partition function. The factors of v can be 
dealt with locally, and the tricky part is to get a handle on the non-local factors Q. To this 
end, it is natural to associate an element of C with each diagram, which will play the role of 
the Boltzmann weight, i.e., the weight of a partially built configuration E' in (2.4). In the 
diagram multiplication d — d 2 • d±, let K,(di, d 2 ) be the number of blocks which are removed 
because they contain no point in the top or bottom rows of d. The non-local part of the 
Boltzmann weight is then Q K ( dl ' d2 ). 

These considerations motivate the definition of the partition algebra [23,36] CAl{Q) as 
the associative algebra over C with basis A L and multiplication defined by 

d 2 d l = Q< d ^{d 2 -d l ). (3.5) 

The partition algebra CA L (Q) can be represented faithfully as an algebra of matrices whose 
rows and columns are indexed by the partition monoid Al. This is indeed the point of 
view that we shall take when constructing the transfer matrix of the flow polynomial and 
manipulating it explicitly (see section 4.2). The elements of Al can then be interpreted as 
the basis states of this representation. With no risk of confusion, we shall therefore use the 
notation A L to refer both to the partition algebra and to the set of basis states. 

We now define a set of generators for the monoid Al- These generators will be the 
elementary building blocks used to define the transfer matrix T l- Apart from the identity 
/, the necessary generators are the join operators Jjj that amalgamate the blocks containing 
points i, j, i' and f, and the detach operators Dj that remove point i' from its block and 
turn it into a singleton. In the pictorial representation this gives rise to the diagrams 

V %' f V V 2' %' V 

-- 1 itHi: i - li i:: n 

1 i j L 1 2 i L 

As a consequence of the above definitions, the product Did within CA L (Q) produces a factor 
Q if the point i' is a singleton in diagram d, and a factor 1 otherwise. In particular, we have 

D 2 t = QDi . (3.6) 
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From these building blocks we can now form the operators representing the addition of 
an edge to the graph that is being built up by the transfer matrix. These are 



Hjj = I + vJij , (3.7a) 
V, = vI + D t . (3.7b) 

The letters H and V stand for horizontal and vertical, where a horizontal edge is understood 
to link vertices within the same layer of the graph (recall Figure 1), and a vertical edge links 
vertices from two adjacent layers. 

Inspecting Figure 1, and labelling the points as in the figure, we can now finally define 
the transfer matrix (with L — k + 1): 

t l = Hoi (n v o H °* ) (n v ( 3 - 8& ) 

\i=k J \i=0 J 

= HoiVoHo2VoHo3---VoH ofc V fc V fc _ 1 ...Vo (3.8b) 
Note the order of indices in the products. 



3.3 Markov trace and eigenvalue amplitudes 

It remains to explain the meaning of the Markov trace Tr in (3.3). The Markov trace 
of any partition diagram d G Al is by definition Q K ^ d \ where n(d) is the number of con- 
nected components in the diagram obtained from d by identifying the points i and i' for 
alH = 1, 2, . . . , L. This identification corresponds to implementing the periodic longitudinal 
boundary conditions in Figure 1. The definition of the Markov trace extends to the partition 
algebra <CAl(Q) by linearity. 

With this definition, (3.3), and (3.8), we are in principle equipped to compute the par- 
tition function Z G ^ n ^(Q,v) as a polynomial in Q and v. A practical problem for going to 
large k is however that the dimension of T L , i.e., the number of basis states A L , grows very 
fast with L: 

dim Tl = \A L \ = B 2L , (3.9) 
where B n are the Bell numbers with exponential generating function (egf) 

00 r> n 

J2^T = exp(e z -l). (3.10) 

n=0 

Considerable progress can nevertheless be made if one takes advantage of the structure 
of the partition algebra [23]. In practical terms this means that the number of points par- 
ticipating in the partitions can be halved from 2L to L. We now explain how this comes 
about. 

Denote by A^ the elements of the partition monoid Al with exactly £ links, and define 
for % — 0, 1, . . . , L 

A L , = \jAf. (3.11) 
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Thanks to (3.4), the Al,% are in fact ideals which, moreover, constitute a filtration of the 
monoid: 

A L ,o C A L ,i C ... C ^ £>£ = /4 L . (3.12) 

This implies immediately that any operator Ol (and T £ in particular) in the partition 
algebra <CAl(Q) has a block-triangular structure with respect to I. The eigenvalues of Ol 
can therefore be found by restricting to . In this restriction, the C- weight of an operator 
that acts on CA^\Q) by diminishing £ will be set to zero. From the point of view of the 
matrix representation of CAl(Q), this restriction amounts to replacing a matrix Ol by 
another matrix 0' L in which all the off-diagonal blocks have been set to zero. Although 0' L 
does not belong to CAl(Q), it is still a well-defined matrix and we can study its eigenvalues, 
which are the same as those of Ol- 

Now £ is conserved, and Ol cannot change the blocks of the partition containing only 
points from the bottom set {1,2,..., L}, since the multiplication has been defined by acting 
on the top points only. Therefore, in the restriction to Af\ L is block-diagonal, with a 
number of identical blocks on the diagonal that corresponds to the number of partitions of 
the bottom points that do not form part of a link. 

As far as the determination of the eigenvalues goes, one can therefore restrict further 
the basis states of A^ to partitions of the top points {1', 2', . . . , 17} only, with precisely £ 
marked blocks (which were the links in the full partition monoid). Note that the marked 
blocks are distinguishable, since the action of Oi, can still exchange their order (relative to 
the now forgotten fixed order of the links with respect to the bottom points {1, 2, ... , L}). 
One can then finally block-diagonalise Ol by rearranging these restricted basis states into 
linear combinations that are irreducible representations A of the symmetric group Se- 

To summarise, all distinct eigenvalues of Ol can be found by studying the irreducible 
representations labelled by £ and A. We have thus the following decomposition of the Markov 
trace 

L 

Tr L = J2J2 ae > x tr ^0L, (3.13) 

e=o \es e 

where now tr^ are ordinary matrix traces. The coefficients o^a are eigenvalue amplitudes, 
which can also be interpreted as the dimensions of the commutant of the partition algebra. 

Consider now A G Sg through its corresponding Young diagram, Y(X) = (X±, A 2 , . . . , A^), 
where Aj is the number of boxes in the i'th row. If there is less than £ rows in Y(X) the 
expression is of course padded with zeros. One then has the result [23, Proposition 3.24] 

^ = ^nw-*-M. (3.i4) 

l - i=0 

We recall that the dimension dim A of the representation A is given by the hook formula [42] 

dimA = ~TTX' (315) 

where h x is the hook length of the box x e Y(X), i.e., the number of boxes to its right, 
plus the number of boxes below it, plus the box itself. We shall sometimes need the total 
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amplitude for a given number of marked blocks. This reads then 

/3 e = ^2 a e ,x dim A. (3.16) 

\es e 

Note that the non-planarity of the graph is crucial for the applicability of (3.13)-(3.14). 
Indeed, had the graph been planar and with periodic transverse boundary conditions (in 
addition to the periodic longitudinal boundary conditions that we assume throughout), the 
transfer matrix could only have changed the cyclic order of the links, and the relevant group 
would not have been Sg, but rather the cyclic group Cg. Its representation theory leads to 
very different expressions [40] for (3.13)-(3.14). A similar remark holds for planar graphs 
with free transverse boundary conditions, in which case the links cannot be permuted at all, 
and only the decomposition according to the value of £ is relevant. 

3.4 Flow polynomial transfer matrix 

The transfer matrix that produces the flow polynomial $G(n,k)(Q) can be taken simply 
as T L of sections 3.2-3.3, i.e., by specialising (3.3) to v — —Q. 

One can however diminish the dimension of the relevant partition algebra by remarking 
that for v = —Q, the vertical operator Vj in (3.7b) is a projector (up to a constant). Indeed, 
by (3.6) one finds V 2 = (vl+ Dj) 2 = v 2 I + (2v + Q)D\, which is a multiple of Vj if and only if 
v = —Q. The normalised projector (— Q) _1 Vj annihilates any partition diagram in which the 
point i' is a singleton. Concerning the reduced partitions of the points {1', 2', . . . , L'} with 
precisely £ marked blocks — as described in section 3.3 — the precise statement is: (— Q)~ l \li 
annihilates any reduced diagram in which %' is an un-marked singleton. 

Since Vj and Vj commute for any the following operator 

L-l 

pl = {-q)- l n v * ( 3 - 17 ) 

i=0 

is also a projector. It annihilates any reduced diagram containing an un-marked singleton. 
We can therefore replace (3.3) by 

Z G{ntk) (Q,-Q) = Tr(T L )"/ fc , (3.18) 

where Tl — PlTl with v = —Q, and consider the trace only over states without un-marked 
singletons. This implies that the flow polynomial [cf. (2.7)] can be finally written as 

*G(n,fc)(Q) = (-l) n g- 2ri Tr(f L )"/ fc , (3.19) 

as the generalised Petersen graph G(n, k) has 3n edges and 2n vertices. The prefactor 
(— l) n Q~ 2n can be absorbed in the definition of the transfer matrix: if we define 

f L = (-l) k Q~ 2k T L , (3.20) 

then (3.19) becomes 

^G(nMQ) = Tr(T L f/ fc . (3.21) 
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The decomposition (3.13)-(3.14) goes through as before, now only with the "no un- 
marked singleton" constraint imposed on the representations labelled by £ and A. 

Remarks. 1. All the entries in the matrix Tj, are polynomials in Q; but this property does 
not hold in general for the matrix elements of Tj,. Some of them may contain terms with 
inverse powers of Q. 5 

2. The structural properties of J l and T £ are obviously the same. 



3.5 Dimensions of representations 

Let us first consider the number of partitions A^ +1 of k + 1 points with £ marked and 
distinguishable blocks. It is given by 

\Af | = k\ [z k ] ((e 2 - if exp (e 2 - 1)) , (3.22) 

as is easily seen by elementary manipulations of the egf of the Bell numbers (the case £ = 0). 
Using (3.10), and the fact that the Stirling numbers of the second kind {^j (or Stirling 
subset numbers) [21] have the following egf [19] 

S\z) = l(e z -lY (3.23a) 

/V| = k\[z k ]S\z) (3.23b) 



we can derive the following closed form for |^4^'' 



'-k 



4"l = «£QHK* (3.24a) 



p=0 

k-i 



s=0 



where we have gone from (3.24a) to (3.24b) by using the well-known expression of the Bell 
numbers in terms of the Stirling subset numbers [19] 



B n = J2\ s \- ( 3 - 25 ) 

It is clear from (3.24) that £\ \ \Af\ = k\, and |^ 0) | = B k . 

Meanwhile, the sum of the eigenvalue amplitudes for a given £ and all possible Young 
diagrams A is given by (3.14)-(3.16): 

i=0 v 7 



5 Let us give an example for k — 3. When we apply the transfer matrix T 4 to the partition (0'1'2'3')(0123), 
we get several partitions with coefficients that are polynomial in Q. In particular, we obtain the partition 
(0'3')(I'2')(0I23) w ith the coefficient (-Q) 5 . If we divide this polynomial by the prefactor (-I) fe Q- 2fe [cf. 
(3.20)], we obtain 1/Q, which is not a polynomial in Q. 
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This is indeed a polynomial in Q; it can be rewritten in terms of the falling factorials 
Q-=IY J= i(Q + ±-j) [21] as: 

P e = E(-WJW- (3.27) 

i=o 

The values we need in this paper are: 

(3o = 1 (3.28a) 

ft = Q-l (3.28b) 

p 2 = Q 2 -3Q + 1 (3.28c) 

p 3 = Q 3 -6Q 2 + 8Q-1 (3.28d) 

/3 4 = g 4 - log 3 + 29g 2 - 24g + 1 (3.28e) 

/3 5 = g 5 - I5g 4 + 75g 3 - 145g 2 + 89g - 1 (3.28f) 

p 6 = Q 6 - 2ig 5 + 160g 4 - 545g 3 + 814g 2 - 415g + 1 (3.28g) 
f3 7 = g 7 -28g 6 + 30ig 5 - 1575g 4 + 4179g 3 - 5243g 2 + 2372g-l (3.28h) 

Just as in the case of the planar partition algebra [see Ref. [31], in particular Eqs. (2.16) 

and (2.20)] the compatibility between the dimensions \A^\ and the amplitudes (3e can be 
expressed in the form of a sumrule: 

This expresses that the number of degrees of freedom per vertex of the graph is indeed Q, 
as expected. We also have the sumrule 

e=o ' r=0 ^ ' 

? (0) 



where the integers B n form the sequence A001861 of [47]. Their egf is 

oo ri(0) n 

= exp (2(e« - 1)) , (3.31) 



ni 

n=0 



as can be deduced from (3.22). 

For fixed k and £, we now introduce the "no un-marked singletons" constraint. We then 
obtain a smaller set of partitions A^ +l of cardinality N k (£) = I^^J. 

The number of partitions of k points with £ marked and distinguishable blocks satisfying 
the "no un-marked singletons" constraint is 

\Af | = k\ [z k ] ((e z - l) e exp (e 2 - 1 - z)) . (3.32) 

The key point is that the number of partitions with no singletons is given by the egf exp(e z — 
1 — z) [19, p. 111]. The numbers associated with this egf are given by 

S n = n![^]exp(e 2 -l-^) = (-lf^T ( H ) B q , (3.33) 
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where the B n are the Bell numbers (3.10). Then, a closed form for the numbers \A^\ reads: 



k 



l^l="EC)g}W (3-34) 



This formula implies that t\ | \A^\, and that \A^\ = k\. 
The sumrule corresponding to (3.29) now reads 



The fact that Q has been replaced by Q — 1 is a manifestation of the "nowhere zero" 
constraint. We also have the sumrule corresponding to (3.30) 



= E^fl = E(*)sf(-1)*-', (3.36) 

£=0 ' r=0 ^ ' 

where the integers are given by (3.30). Their egf is 

00 5(0) n 

J2^T- = exp(2(e 2 -l)-z) . (3.37) 



n=0 



Notice that the difference \ A^ \ — \A^ \ gives the number of partitions of the set {1, 2, . . . , k} 
with £ marked points and with at least one un-marked singleton. These partitions do not 
contribute to the final result, as they are associated to null eigenvalues. 



4 Flow polynomial for the generalised Petersen graphs 

4.1 General theory 

Let us start with the simplest case k — 1. The graph G(n, 1) is isomorphic to a cyclic lad- 
der of width 2. From the known Potts-model partition function [18, and references therein], 
one can easily derive 

*G(n,i)(Q) = (Q 2 -3Q + l)(-ir + (Q-l)(Q-3r + (Q-2r, (4.1) 

where the eigenvalues /j, — —1,Q — 3, Q — 2 correspond to the sectors with £ = 2, 1, links, 
respectively. 

Let us now focus on k > 2. The flow polynomial (3.21) for the generalised Petersen graph 
G(nk, k) can be written as 

fc+i 

®G(nk,k)(Q) = $^$^,Atr^ A (T fc+1 ) n , (4.2) 

e=o xes e 
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where the amplitudes are given by (3.14). This formula is the most general one. In 
terms of the non-zero eigenvalues Hk,e,\,s of the transfer matrix T^+i, it reads: 

fc+l N k (i) 

®G(nk,k)(Q) = ^ J2 Vll,\,s , ( 4 -3) 

e=o xeSi s=i 

where N k (£) = (l^jj [cf. (3.32)/(3.34)]. 

The flow polynomial $G(nfc,fc) (4-3) is obtained as a linear combination of ordinary matrix 
traces (i.e., sums of powers of the eigenvalues /ik,e,\,s) with definite coefficients (3.14). 
This is all what we need to compute <&G{nk,k) rigorously from the various diagonal blocks 
of Tx,. It is worth stressing that Eq. (4.3) holds true irrespective of whether some eigen- 
values happen to be identical or not. However, it will turn out useful in practice to have a 
"complete" decomposition of the flow polynomial, in the sense that there are no equalities 
among eigenvalues (except for what is accounted for in the multiplicities eee,\)- This reason 
motivates a deeper study of the eigenvalue structure of the transfer matrix T L . 

To simplify the notation, we will denote T k+1:i the diagonal block of the full transfer 
matrix T fc+1 corresponding to partitions with exactly £ links. We will denote Tk+i,e,\ the 
diagonal block of Tk+i/ corresponding to the irreducible representation A of the group Si. 
Similar notation will be used for the diagonals blocks of the transfer matrix T^+i- 

For 1 < k < 7, our exact computations (described in detail in section 4.2) have revealed 
the following structure described in Lemmas 4.1, and 4.2, Corollary 4.3, and Lemma 4.4 
below: 

Lemma 4.1 Fix 1 < k < 7. Then for any < i < k + 1, the diagonal block Tk+i,e. contains 
the trivial eigenvalue Hk,k+i — ( — l) fe with multiplicity A r fc l (£)£! 7 where 

and S n is given in (3.33). 



Remarks. 1. We have proven that this lemma can be extended to any k > 1. This general 
result corresponds to Lemma A.l in Appendix A. 

2. The eigenvalue Hk,k+i = {— l) fe will be called "trivial" in the following, because it 
is directly related to the structural properties of the transfer matrix J l- In the proof of 
Lemma A.l, we show that, when we order the basis states appropriately, the diagonal block 
Jk+i,£ has a block triangular structure with two diagonal blocks: 

• One block is a diagonal matrix of dimensionality Nk,i(£)£\ with all diagonal entries 
equal to the trivial eigenvalue fj,k,k+i- 

• The other block is not in general a diagonal matrix, and gives rise to the "non-trivial" 
eigenvalues. 
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3. In particular, the general Lemma A.l imply that for any k > 1 and £ = k + 1, the 
multiplicity of the trivial eigenvalue fi k k+i is given by N k (k + 1) = N k i{k + 1) (k + 1)! = 
(k + l)\. 



Lemma 4.2 Fix 1 < k < 7, and £ = k. Then, there are k non-trivial eigenvalues fi k ,k,s 
in the diagonal block Tk+i,k (once we have removed the trivial eigenvalue fj,k,k+i = ( — l) k )- 
Each of these eigenvalues fj,k,k,s has multiplicity k\. 



Remarks. 1. We have proven that this lemma can be extended to any k > 1. This general 
result corresponds to Lemma A. 2 in Appendix A. 

2. For each irreducible representation A of the group S k , the diagonal block T k+ i tkt \ 
contains the same k non-trivial eigenvalues /ik,k,s- 



Corollary 4.3 Fix 1 < k < 7. Then for all pairs (£, A) with < £ < k — 1 and A G Se, 
the diagonal block Tk+i,t,\ has N kj0 (£) non-trivial eigenvalues [ik,e,\,s (once we have removed 
the trivial eigenvalue fik+i,k = (~ which appears for every 1 < £ < k + 1). The value of 
Nk,o(£) is given by: 



Sk+i . ) s k - P 

[Sk-p + Sk+i-p] 



k 

p=l 

k 





p=0 

k\ \p 

V 



if£ = 

if\<£<k-\ 

if£ = k 
if£ = k + l 



(4.5) 



where S n is given by (3.33), and Y e is the number of Young tableaux with £ cells. 



Remarks. 1. The number Ye of Young tableaux with £ cells (sequence A000085 in [47]) 
has the following egf: 

Y e = ^dimA = £\ [z e ] exp (z + z 2 ) (4.6a) 

xes e 

» - §SC) 

2. This corollary follows from Lemmas 4.1 and 4.2. As these two lemmas can be gener- 
alized for any k > 1, the more general Corollary A. 3 in Appendix A follows. 



Lemma 4.4 Fix 1 < k < 7. Then for each pair (£, A) with < £ < k and A G Si, the 
non-trivial eigenvalues Hk,e,\,s ( or A^.M when £ = k) are all distinct and also different from 
the trivial eigenvalue Hk,k+i — (~^) h ■ 
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Remarks. 1. We do not know how to prove the extension of this result to larger values of 
k > 8; but we conjecture that Lemma 4.4 holds true for any k > 1. 

2. The fact that all eigenvalues Hk,t,\,s are distinct in the generic case, implies that we do 
not have to worry about the non-degenerate-dominance condition in the Beraha-Kahane- 
Weiss theorem (see section 5). 

Lemmas 4.1, 4.2, and 4.4 show that for 1 < k < 7, we have found that the decomposition 
(4.3) is " almost complete" , as all eigenvalues /ik,e,x,s are distinct, except for two cases: 

(a) The trivial eigenvalue fXk,k+i — (— l) fc > that appears in every sector with 1 < t < k + 1. 

(b) The eigenvalues /ik,k,s = Hk,k,\,s (°f the sector t = k) do not depend on the value of the 
representation A G Sk- 

This means that, if we want to have a "complete" decomposition of the flow polynomial, 
then the amplitudes for the eigenvalues fik,k+i and Hk,k,s are n °t longer given by a single co- 
efficient a^; but by a linear combination of several coefficients. In particular, the coefficient 
associated to the eigenvalues f/,k,k,s should be equal to (3k (3.16)/(3.26). The coefficient 7^+1 
associated to the trivial eigenvalue /ik,k+i should be given by 

k 

lk+ i = Pk+i + Y,Nk,i(e)Pe, (4.7) 

where Nk,i(£) is given in (4.4). 

Lemmas 4.1, 4.2, and 4.4, Corollary 4.3, and the above discussion imply that the flow 
polynomial for the generalised Petersen graph G(nk, k) is given by the following "complete" 
decomposition: 

Corollary 4.5 Fix 1 < k < 7. Then the flow polynomial for the generalised Petersen graph 
G(nk, k) is given for any n>\ by 

k-l iVfc.oW k 

^nkMQ) = EE Q ^ E <^+^E<M+7 fc+ i(-i) nfc (4-8) 

e=o xes e s=i s=i 

where is given in (3.26) / '(3.27) , 7^+1 is given by (4-7), and all eigenvalues {^k,£,x,s, Hk,k,s, 
(— l) fc } are distinct. 

Remark. We conjecture that the above "complete" decomposition for the flow polynomial 
$G(nfc,fc) holds true for any k > 8. 

We can compute the coefficients 7^+1 by using (4.7), the expressions for the amplitudes 
(3 e (3.28), and the values of N k M (4-4). The results for 1 < k < 7 are: 



72 


= Q 2 


- 3Q 2 + 1 






(4.9a) 


73 


= Q 3 


- 5Q 2 + 6Q-1 






(4.9b) 


74 


= Q A 


- 7Q 3 + 15Q 2 - 


ng + i 




(4.9c) 


75 


= Q 5 


- 9Q 4 + 28Q 3 - 


38Q 2 + 20Q - 1 




(4.9d) 


76 


= Q 6 


-ng 5 + 45g 4 - 


- 90g 3 + 90Q 2 - 


27Q+ 1 


(4.9e) 


77 


= Q 7 


- 13Q 6 + QQQ 5 - 


- 175Q 4 + 260Q 3 


- 207Q 2 + 70Q - 1 


(4.9f) 


78 


= Q s 


- 15Q 7 + 91Q 6 - 


- 30ig 5 + 595Q 4 


- 707Q 3 + 469Q 2 - 135Q + 1 


(4.9g) 
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For instance, for k = 3, we have that 74 = 0i + 4/3 2 + 3/3 3 + /3 4 . 

The total number of distinct eigenvalues as a function of k is given by 



k 



(4.10) 



where the N kfi (£) are given in (4.5). The values are D k = 3,7,36,229,1658,12803,105934 
for k = 1,...,7. 

Remarks. I. It is important to notice that, once the eigenvalues of the transfer matrix 
T k+ i [or those of T k+ i (3.20)] are obtained, as well as their multiplicities, Corollary 4.5 
establishes the exact form of the flow polynomial <&G(nk,k) for any 1 < k < 7 and any n > 1. 
In particular, the exactness of the analysis does not depend on the validity of Lemma 4.4 
for larger values of k > 8. 

4.2 Practical procedure 

We have written a perl script to compute the symbolic transfer matrix T^+i using ideas 
similar to those already explained in [26,27]. For 1 < k < 4, we have checked our programs 
using mathematica. Further checks were performed with code written in C that allows us 
to compute the leading eigenvalue for given values of k, £, and A = (£) — the completely 
symmetric irreducible representation of Se- 

The first step is to obtain the relevant diagonal blocks T k+ i^ \ of the transfer matrix 
J k+ i. We first fix the value oi < £ < k + 1, and a bottom- row configuration compatible 
with the chosen value of £ and the "no un-marked singletons" condition. We then determine 
the basis of the relevant partition space of dimensionality N k (£) [cf. (3.34)]. Indeed, the 
result does not depend on the chosen bottom- row partition. 6 

We now choose an irreducible representation A^of the symmetric group Se of dimensional- 
ity dim A. To obtain the relevant diagonal block T k+ ± t e,x corresponding to A, we simply take 
as our basis vector those linear combinations of the "standard" basis with the appropriate 
properties under Se- Finally, in order to "extract" the trivial eigenvalues 7 Q 2k for £ > 1, we 
exploit the structure of T k +i,t,\ found in the proof of Lemma A.l. 

In this way, we obtain all the diagonal blocks T k+1 t e,\ with 1 < A; < 7, < £ < k + 1, and 
all the A G Se- We then have to compute the traces 



for integer values of n > 1. For small values of k, this computation can be made symbolically 
for any n using a symbolic algebraic manipulator program, such as mathematica. However, 

6 To be more precise, there are bottom-row configurations that give rise to a basis of dimensionality pN k {£) 
[cf. (3.34)] for some integer p > 2. In these cases, one indeed obtains the same N k (£) eigenvalues; but each 
one repeated p times. 

7 Plcase note that the trivial eigenvalue for the matrix T^+i (3.20) is given by /i k ,k+i = (~ l) fe - But the 
corresponding eigenvalue for the matrix Tfc + i is (— l) fc (— l) k Q 2k = Q 2k . 



N k , (e) 




(4.11) 



s=l 
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for larger values of k (say, k — 6, 7) this is not feasible, as we have blocks of dimension as 
large as 11 816 (for k = 7, £ = 3, and A = (2, 1)). 

A key issue in the subsequent analysis is the dependence of the above traces (4.11) on 
Q. The needed information is given by the following Lemma (the proof can be found in 
Appendix A): 

Lemma 4.6 Let k>l,0<£<k + l, and A be an irreducible representation of Sg. Then 
trT£ +lj£)A [cf. (4-11)] is a polynomial in Q of degree at most n[k + min(l — £, 0)]. 

We have made use of the following tricks: 

• We know that the traces (4.11) are polynomials in Q of degree d = n[k + min(l — £, 0)] 
by Lemma 4.6. Therefore, we only have to compute the evaluation of each trace at 
d+1 integer values of Q ^ 0, and then reconstruct the corresponding polynomial using 
Lagrange's interpolation method. In order to check the result, we always compute at 
least d + 2 values of the trace. Please note that we compute the evaluation of the trace 
(4.11) by first computing the evaluation of the trace trT£ +1 ^ A , and then multiplying 
the result by the factor (-l) nk Q^ 2nk [cf. (3.20)]. Indeed, this factor is not defined for 
Q = 0. 

• Not all the entries of the transfer matrix Tk+i,e,x are integers for integer values of Q; 
rather they are rational numbers. As we want to perform the trace computation with 
(infinite-precision) integer arithmetic, for each value of Q, we multiplied the matrix 
by the minimum (positive) integer value such that all entries are integers. After the 
computation is done, we reconstructed the true solution by dividing by the appropriate 
factor. 

• The integers involved in the actual calculations are very large. Therefore, we compute 
the value of the trace for a given value of Q ^ using modular arithmetic for a 
given set of prime numbers p < 65521 < 2 16 (we need up to 65 different primes). We 
then reconstruct the value using the Chinese remainder theorem using infinite-precision 
arithmetic in mathematica. We always use at least one more prime than needed, in 
order to check the result. To accelerate the computation of the trace using modular 
arithmetic, we use a program written in C. 

• For 1 < k < 6, we are able to compute the traces for many different values of n. 
However, for k = 7, the computation is so demanding, that we have focused on powers 
of the type n = 2 q + 1 with integer q > 1. The reason why we consider odd powers for 
k = 7 will become clear in section 6. 

Once the traces are computed, we can form the flow polynomial using (4.8). Notice that 
for 1 < k < 7, everything in this formula is exactly known. The zeros of the flow polynomials 
are obtained using the program mpsolve [12,13]. This software has the advantage that if 
one requests the zeros with 50 digit precision (as in our case), the results are guaranteed to 
have at least such precision. 

Remark. To give a clear idea of what has been achieved, consider the case k = 7, and 
more specifically the computation of $ G for the graph G = G(17k,k), which is the largest 
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computation undertaken in this work. Applying naively (3.8) within the diagrammatic basis 
would imply computing the 17th power of T^ +1 , which according to (3.9) is a matrix of 
dimension B 2 (k+i) — 10 480 142 147 whose entries are polynomials in Q of degree 2k. The 
decomposition of T fc+1 and use of the "no un-marked singleton" constraint has reduced the 
computation to the sum over 31 blocks, the largest of which has dimension 11 816. Even with 
these tricks, the computation took around six months, using 50-80 processors, corresponding 
to some 30 years of CPU time. 

4.3 Additional checks 

Because our results are derived using software, we have performed some tests in order 
to ensure that the results are right. First of all, for the smallest members of each family 
G(nk,k), we have computed the flow polynomial using maple. Indeed, both results agree 
in all cases. 

For the cubic graphs G(n, k) that we are considering, we may improve on (2.12) by adding 
a few more terms: 

aWQ) = Q lEHVl+1 - \e\qW-w + ( ™ E 2 I-1) - Q'*'-'^- 1 

)(|E| ~ 2) - \V\(\E\ - 2)^ Q\ E \-\ y \~ 2 + .... (4.12) 

In this expression, the coefficient of Q\ E \~\ V \ _1 arises from two contributions in which E\E' is 
respectively two edges, and three edges all incident on the same vertex. The contributions to 
the coefficient of Q\ E \~\ V \~ 2 are slightly more complicated to characterise. Inserting |V| = 2n 
and \E\ = 3n we obtain 

*G(n,*)(Q) = Q n+1 (l " SnQ- 1 + |(9n - 7)Q- 2 - ^(3n - 2)(3n - 5)Q~ 3 + ...). (4.13) 

We have checked that for all the graphs G(nk, k) we have considered, the flow polynomials 
obtained from the procedure outlined above indeed satisfy (4.13). 

There are some theorems that give us some information about the location of the real 
zeros of the flow polynomial. The first theorem applies to a general bridgeless graph; while 
the second one applies only to cubic graphs (i.e., it is valid for G(nk, k) with n > 2): 

Theorem 4.7 (Wakelin [55]) Let G be a bridgeless graph with \V\ vertices, \E\ edges, b 
blocks, and no isolated vertices. Then: 

• &g(Q) is non-zero with sign (— l)^-^!-*- 1 for Q G (— oo, 1); 

• $g(Q) has a zero of multiplicity b at Q — 1. 

• $g(<5) is non-zero with sign (-l)l s H y l+ 6 + 1 forQ G (1, p]. 

Theorem 4.8 (Jackson [24, 25]) Let G be a 3-connected cubic graph with \V\ vertices and 
\E\ edges. Then: 

• $g(<5) is non-zero with sign (-1)1-^1-^1 for Q G (1,2); 
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• &g(Q) has a zero of multiplicity 1 at Q = 2. 

• $g(Q) is non-zero with sign (— l^H^I+i for Q E (2,5), where 5 ~ 2.546 is the flow 
root of the cube in the interval (2, 3) [i.e., the zero in this interval ofQ 3 — 9Q 2 +29Q— 32/. 

The generalised Petersen graphs G(nk,k) satisfy |V| = 2nk, \E\ = 3nk, and 6 = 1. 
Therefore, the above theorems imply that, for any n > 2, has only two simple real 

roots in the interval (-oo,<J): Q = 1,2. The sign for Q e (-oo, 1)U(2,«J) is that of (-l) nfc+1 , 
and it has the opposite sign for Q e (1,2). For large enough Q > 0, the sign of $g(Q) is 
always positive. Therefore, for even k and any n > 3, or for odd k and even n > 4, this 
implies the existence of a real zero in [5, oo). 8 

We have explicitly checked that all the computed flow polynomials $G{nk,k) have only 
two simple roots in the whole interval (—oo,S), namely Q — 1,2. Furthermore, for all even 
(resp. odd) k, and all n > 3 (resp. all even n > 4), the polynomial $G(nfc,fc) has at least one 
root in [S, oo). 

If G is a cubic graph, one can easily see whether Q = 3 is a flow root or not [17, 
Proposition 6.4.2]: 

Theorem 4.9 ^4 bridgeless cubic graph G has a nowhere zero 3-flow if and only if it is 
bipartite. 

As the graphs G(nk, k) are bipartite if and only if k is odd and n is even, Theorem 4.9 
implies that $G(nk,k)(Q) has at least one factor Q — 3 whenever G(nk,k) is not bipartite: 
namely, when k is even or n is odd (or both). When G(nk, k) is bipartite, then $G(nfc,fc)(3) > 
0. We have checked these facts in all the flow polynomials we have explicitly computed in 
this work. 

5 The Beraha— Kahane— Weiss theorem 

A central role in the subsequent analysis (see sections 6-7) is played by a theorem on 
analytic functions due to Beraha, Kahane and Weiss (BKW) [8-11] and generalised slightly 
by Sokal [48]. The situation is as follows: Let D be a domain (connected open set) in the 
complex plane, and let a±, . . . , %, /ii, . . . , /im (M > 2) be analytic functions on D, none of 
which is identically zero. For each integer n > 0, define 

M 

f n (z) = Y, a k{z) ^ k {z) n . (5.1) 
k=\ 

We are interested in the zero sets 

Z(fn) = {zeD: f n (z) = 0} (5.2) 



8 The lower bound is sharp, as the cube is isomorphic to the generalised Petersen graph G(4,3), which 
has a zero at Q = 5. 
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and in particular in their limit sets as n — > oo: 

liminf Z(f n ) = {z G D: every neighbourhood U 3 z has a nonempty intersection 

with all but finitely many of the sets Z(f n )} (5.3) 

limsup Z(/ n ) = {z G D : every neighbourhood U 3 z has a nonempty intersection 

with infinitely many of the sets Z(f n )} (5.4) 

Let us call an index k dominant at z if \fi k (z)\ > \fJ>i(z)\ for alH (1 < I < M); and let us 
write 

D k = {z G D: k is dominant at z} . (5.5) 
Then the limiting zero sets can be completely characterised as follows: 

Theorem 5.1 (Beraha— Kahane— Weiss [8—11,48]) Let D be a domain in C, and let 
oti, . . . ,otM, fJ>i, ■ ■ ■ , I^m (M > 2) be analytic functions on D, none of which is identically 
zero. Let us further assume a "no -degenerate- dominance" condition: there do not exist in- 
dices k 7^ k' such that fj, k = cofi k ' f or some constant oj with \oj\ = 1 and such that D k (— D k i ) 
has nonempty interior. For each integer n > 0, define f n by 

M 

fn(z) = ^2a k (z) n k {z) n . 
k=l 

Then liminf Z(f n ) = limsupZ(/ n ) 7 and a point z lies in this set if and only if either 

(a) There is a unique dominant index k at z, and a k (z) = 0; or 

(b) There are two or more dominant indices at z. 

Note that case (a) consists of isolated points in D, while case (b) consists of curves (plus 
possibly isolated points where all the \i k vanish simultaneously). Henceforth we shall denote 
by B the locus of points satisfying condition (b). 

We shall often refer to the functions fi k as "eigenvalues" , and to the a k as "amplitudes" , 
because that is exactly how they arise in the transfer matrix formalism. 

In Ref. [11, p. 55], Beraha, Kahane, and Weiss give (without proof) the following corollary 
based on their proof of Theorem 5.1: 

Corollary 5.2 If z is a real isolated limiting point, the polynomials f n have real coefficients, 
the dominant eigenvalue /x* and its coefficient a+ are both real in an interval (z — e, z + e) 
for some e > 0, and if a[(z) ^ 0, then z is the limit of a real sequence {z n } with f n (z n ) = 0. 

The proof is simple: For any sufficiently small e > 0, a±(z ± e) are non-zero of opposite 
sign, and still dominant. Then, for all sufficiently large n (depending on e), f n (z ± e) are 
real of opposite sign. Therefore there exists a root in between. 

In the next section we need some general results on the existence of a real sequence of 
zeros {z n } such that it converges to a real non-isolated limiting point. These results can be 
summarised in the following 
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Lemma 5.3 Let us suppose that z is a real non-isolated limiting point, such that exactly two 
dominant eigenvalues jii and /i 2 become equimodular at z. Let us further suppose that: 

(a) The two dominant eigenvalues are analytic functions in a neighbourhood of z. 

(b) The eigenvalue Hi (resp. fi 2 ) is dominant (resp. subdominant) in the interval [z — e,z), 
and is subdominant (resp. dominant) in the interval (z, z + e] for some e > 0. 

(c) The corresponding amplitudes a\ and a 2 do not vanish at z. 

(d) The eigenvalues fii,fi 2 and the amplitudes a\,a 2 are real in a real neighbourhood of z. 
Then, 

1. If a\(z)a 2 (z) > and fii(z) = —ji 2 (z), then for oddn, there is a real sequence of zeros 
{z n } converging to z. 

2. If ai(z) a 2 (z) < and fii(z) = —fi 2 (z), then for even n, there is a real sequence of 
zeros {z n } converging to z. 

3. If ai(z)a 2 (z) < and fii(z) = fi 2 (z), then for all n, there is a real sequence of zeros 
{z n } converging to z. 



Proof. The function /„ (5.1) can be written as 



N 

i=3 



(5.6) 



where the N — 2 other eigenvalues are subdominant in a neighbourhood of z. Then, for every 
sufficiently small e > we can find real numbers < n, r 2 < 1 such that 



for all i > 2 and w e [z — e,z), and 



Hi(w) 



< n, 



(5.7) 



fi 2 (w) 



< r 2 , 



(5.1 



for alH = 1 or i > 3 and w G (z, z + e] . Let M the maximum of the amplitudes cti in the 
interval [z — e,z + e]. We choose e small enough so that the signs of the dominant eigenvalues 
and amplitudes is the same as at z. Then, we have that 



f n (z-e) = nl 
f n (z + e) = & 



N 



a i 



(\ n N / \ n 



(5.9a) 



(5.9b) 



Then, for large enough n, the quantities f n (z ± e) have opposite sign in the following cases: 
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1. a.\OLi > and (ni/ii 2 ) n < 0, which occurs if /ii/i 2 < and n is odd. 

2. 0i0 2 < and (fi 1 /fi 2 ) n > 0, which occurs if iii/i 2 > (then n can have either parity) 
or if if /Zi/i2 < and n is even. 

In any of these cases, the continuous function /„ attains values of distinct signs at the 
endpoints of the interval [z — e, z + e] , therefore there should be a zero at some point inside 
this interval. I 

Remarks. 1. In the fourth case 0:10:2 > and ^ = /i 2 , the zeros converging to z are 
non-real. 

2. If the derivative of the ratio fi\/fi 2 is nonvanishing at z, then condition (b) [or the 
same condition with ^ and /i 2 interchanged] necessarily holds. However, the converse is 
false: it is possible for condition (b) to hold even if fi\/fi 2 has a vanishing derivative at z. 

As remarked on Ref. [11, p. 55], the convergence ratio for isolated and non-isolated 
limiting points is rather different: exponentially fast for the former \z n — z I < Ar n , and 
\z — z n I < An' 1 for the latter, as n —> 00. 

6 Real zeros of the flow polynomials <&G(nk,k) 

In this section we will discuss the real zeros of the flow polynomials &G(nk,k) for 1 < A; < 7. 
In particular, we will focus on the real zeros around Q — 5, and on the existence of real zeros 
Q > 5. 

6.1 k < 5 

For 1 < k < 3, we compute the flow polynomials of G(nk, k) and their roots for all n 
in the range 1 < n < 30. All the flow roots we have found are smaller than Q = 4. We 
conjecture that this holds for larger n as well. For k = 4, 5 we find flow polynomials with 
real roots greater than Q = 4: 

• G(28,4) has two real roots greater than Q = 4: Q x w 4.0002086861 and Q 2 w 
4.3876416603. As n grows, the maximal real root of $G(4n,4) tends to the value 
Q c (4) w 4.5697435537. 

• G(30,5) has two real roots greater than Q = 4: Q x w 4.0000786673 and Q 2 ~ 
4.4867394006. As n grows, the maximal real root of $G(5n,5) tends to the value 
Q c (5) w 4.9029018077. 

Both families thus provide counter-examples to Welsh's conjecture (cf. Conjecture 1.3). 

Remark. The values of Q c {k) mentioned above are obtained by determining the eigen- 
value crossing that corresponds to case (b) of Theorem 5.1. See Conjecture 7.2 below for a 
more precise statement about the leading eigenvalues at Q c (k). 
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6.2 k = 6, 7: The isolated limiting point Q = 5 



The family G(6n, 6) also provides a counter-example to Conjecture 1.3: it displays a 
real zero converging to Q = 5 from below; for instance, G(144, 6) has a real zero at Q « 
4.9987003379. 

This empirical observation will be made rigorous by applying Corollary 5.2 to the family 
G(6n, 6) at Q — 5. In this case, there is a unique (and therefore, real) leading eigenvalue, 
which corresponds to £ = 3 and A = (3): /i 3 ,( 3 ) iVc (5) « 177.122. The corresponding amplitude 
«3,(3) is given by (3.14), 

"3,(3) (Q) = ^Q(Q-l)(Q-5), (6.1) 

and is a polynomial in Q with a single zero at Q = 5. Therefore, the BKW theorem implies 
that Q = 5 is an isolated limiting point for this family. The fact that the limiting point 
Q = 5 is isolated implies that there is a interval of radius e around Q = 5 where there are 
no other limiting points. Therefore, in this interval the eigenvalue ^3,(3),* is still dominant, 
and hence real. Therefore, Corollary 5.2 implies that there is a sequence of real zeros {Q n } 
converging to Q = 5. We can go a little bit further and show that: 

Corollary 6.1 The point Q = 5 is an isolated limiting point for the family G(6n, 6) with 
n > 3. There is a sequence of real zeros {Q n } of the flow polynomial $G(6n,6) that converges 
to Q = 5 from below. 



PROOF. Let us consider the point Q = 5— e with e small enough so that: 1) the leading eigen- 
value A*3,(3),*(Q) > 0; 2) all other N sub-leading eigenvalues /ij satisfy |/^(Q)//i 3 ,(3),*((3)| < 
r < 1; 3) the sub-leading amplitudes |<x,| < M are bounded; and 4) 0:3,(3) (Q) < 0. We 
can always choose e so that these conditions are fulfilled, as the eigenvalues and ampli- 
tudes are analytic functions of Q in a neighbourhood of Q = 5. In addition, we know that 
$G(6n,6)(5) > 0, due to Lemma 3.1. Therefore, we only need to show that $G(6n,6)(Q) < 
for n large enough. This is easy as, 



$G(6n,6)(Q) = A*3,(3),*(Q)' 



o 3)(3) (g) + Voj(g) 



TV 



A*3,(3),*(Q) 



(6.2) 



We can always choose iVo such that NMr N ° < |o 3 ,( 3 )(5 — e)|. Then, for all n > N , the 
sign of $G(6n,6)(Q) is that of 03,(3) (Q) (i- e -> negative), so that there should be a point Q in 
(5 - e, 5) such that ® G (6n,6)(Qo) =0. ■ 

The same situation applies to the family G(7n, 7) at Q = 5. In this case, the unique (and 
real) leading eigenvalue is // 3 ,( 3 ),*(5) ~ —621.779. Therefore, Q = 5 is an isolated limiting 
point for this family, and there is a real sequence of zeros {Q n } converging to Q — 5. 
Because /^ 3 ,( 3 ),*(5) < 0, the above arguments imply that the convergence of the sequence 
{Qn} to Q = 5 is a bit more complicated: 

Corollary 6.2 TTie pom£ Q = 5 is an isolated limiting point for the family G(7n, 7) wrat/j 
n > 3. There is a sequence of real zeros {Q n } of the flow polynomial $G(7n,7) converging to 
Q = 5. TTie sub-sequence with odd n (resp. even n) converges to Q = 5 /rom ak>ve (Vesp. 
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Proof. If we consider even n, then the proof is as before; therefore the sequence {Q2 P } P eN 
converges to Q = 5 from below. On the contrary, for odd n, then the sign of the leading 
term tt3 i (3)/i3 i (3) i *(<5) n is positive for Q < 5, and negative for Q > 5. Therefore, a trivial 
modification of the above arguments leads to the convergence of the sub-sequence {Q2 P +i} p gn 
to Q — 5 from above. I 

Indeed, Corollaries 6.1 and 6.2 imply Theorem 1.7. Moreover, Corollary 6.2 implies the 
existence of real flow roots arbitrarily close to Q = 5 on both sides (above and below). 

6.3 k = 7: Real flow— polynomial zeros larger than Q = 5 

There are also real non-isolated limiting points for the families G(nk, k) with 1 < k < 
7. These non-isolated limiting points corresponds to "crossings" between two dominant 
eigenvalues, each of them coming from a block of the transfer matrix Tk+i,t,\ with different 
values of £. Therefore, both dominant eigenvalues are real and analytic in some real interval 
around the non-isolated limiting point Q c (k). In some cases, there corresponds a sequence 
of real zeros {Q n } converging to that limiting point. 

The family G(6n, 6) would be in principle a good candidate for having real roots larger 
than Q = 5: a direct calculation shows that Q C (Q) ~ 5.1079785012 is a non-isolated limiting 
point for this family. However, the actual computation of all the members of this family up 
to G(144, 6) does not reveal any zero Q > 5. 

The explanation is simple: at Q = Q C {Q) both dominant eigenvalues are equal Ai 5 (i)^((5c(6)) 
— ^3,(3),*( ( 5c(6)) ~ 169.757, and the corresponding amplitudes are both positive. Therefore, 
the hypotheses of Lemma 5.3 are not satisfied, and we cannot find real zeros converging to 
the non-isolated limiting point Q c (6). They should be non-real. 

Finally, we have found that the family G(7n, 7) with odd n has members with the desired 
property. Even though G(63, 7) does not have any zeros larger than Q = 5, the next member 
that we have computed (2(119, 7) has two zeros: Qi w 5.0000197675, and Q 2 « 5.1653424423. 
The flow polynomial for (2(119, 7) has degree 120 and can be written as 

*G(ii9,7)(Q) = (Q-1)(Q-2)(Q-3)P 117 (Q), (6.3) 

where P U7 (Q) = Q 117 - 351Q 116 + 61191Q 115 - 7064107Q 114 + ... is a polynomial in Q of 
degree 117. The coefficients of this polynomial are given in Table 2. We can formalise the 
existence of such real roots greater than Q = 5 in the following way: 

Proposition 6.3 Let $g(Q) be the flow polynomial of the generalised Petersen graph G = 
G(119, 7). Then it has a real zero in the interval (5 + 10 -5 , 5 + 2 x 10~ 5 ) ; and another real 
zero in the interval (516534 x 10" 5 , 516535 x 10" 5 ). 

Proof. We can evaluate the polynomial <&q(Q) at the two end-points of the interval (5 + 
10~ 5 , 5 + 2 x 10~ 5 ) using exact rational arithmetic and find results of distinct sign: $g(5 + 
10~ 5 ) w +2.21791 x 10 42 , and $ G (5+2x 10~ 5 ) w -5.27937x 10 40 . Therefore, the intermediate 
value theorem ensures the existence of a zero of $ G m the open interval (5 + 10~ 5 , 5+2 x 10~ 5 ). 

The same procedure can be carried out for the second interval: $g (516534 x 10~ 5 ) w 
-1.46592 x 10 42 , and $ G (516535 x 10~ 5 ) w +4.53729 x 10 42 . I 
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Remark. As a curiosity, $ G (ii9,7)(5) = 4488918995790513676672232799446257724715600, 
and $ G (ii9,7)(4) = 1133172760943853528. 

As a matter of fact, we also expect that the family G(7n, 7) has, for every large enough 
odd n, two real zeros larger than Q = 5: one converging to Q = 5 from above (as we 
have discussed in the previous section), and the other one converging to the limiting point 
Q c (7) ~ 5.2352605291 from below. The complete statement is given by the following 

Lemma 6.4 The point Q c (7) G (5235260 x 10 -6 , 5235261 x 10~ 6 ) is a non-isolated limiting 
point for the family G{7n, 7) with n > 3. There is a sequence of real zeros {Q n } of the flow 
polynomial $G(7n,7) converging to Q c (7) from below for odd n. 

Proof. In the interval [5235260 x 10~ 6 , 5235261 x 10~ 6 ], we only find two dominant eigen- 
values Ao,* and A 3i ( 3 ) ; *, each of them coming from a different £ sector. Therefore, both 
eigenvalues and its sum A 0i * + A 3i ( 3 )^ are analytic functions of Q in this interval. If we 
evaluate A ,* + A 3) ( 3 )^ at the two end-points of this interval using mathematica, we find 
results of distinct sign: A ,*(5235260 x 10~ 6 ) + A 3)(3)i * (5235260 x 10~ 6 ) « 0.000917, and 
A ,*(5235261 x 10~ 6 ) + A 3)(3) ^(5235261 x 10~ 6 ) « -0.000817. Therefore, there exists an 
intermediate value Q c {7) in the open interval (5235260 x 10~ 6 , 5235261 x 10 -6 ) such that 
Ao,*(Qc(7)) + A 3i ( 3 ) iVr (Q c (7)) = 0. Exactly at Q = Q c (7), these two eigenvalues have opposite 
signs A ,*(Q C (7)) = — A 3i ( 3 )^(<5c(7)) ~ —565.833. The corresponding amplitudes are both 
positive for any Q > 5. Therefore, according to Lemma 5.3 we find real zeros converging to 
Qc{7) only for large enough odd values of n. 

The fact that the convergence if from below comes from the fact that ao^(Q c (7)) — 
«3,(3),*(Qc(7)) > 0, so the sign of the flow polynomial at Q c (7) is positive $G(7n,7)(Qc(7)) ~ 
K*(<2c(7)) - a 3i(3) ^(Q c (7))]A ^(Q c (7)) ri > 0. However, its sign at Q c (7) - e (for small 
enough values of e > 0) $ G (7n,7)(Q c (7) - e) « a 3i(3) ^(Q c (7))X 3i(3) ^(Q c (7)) n < for odd n. 
Therefore, there should be a root in between. | 

Proposition 6.3 and Lemma 6.4 imply Theorem 1.8. 

7 Numerical results 

We now turn to the computation of the limiting curves For any fixed integer k > 1, 
the curve is defined as the set of non-isolated limiting points of the zeros of the flow 
polynomial $G(nk,k), as n — > oo. We have computed these curves for k = 1,2, ... ,7 using 
the direct-search method, as described in [26]. The results are shown in Figures 2-5, both 
in the complex Q-plane and in the complex l/Q-plane. 

We have noticed several empirical patterns in Figures 2-5. The first one is related to the 
number of outward branches of the limiting curve and the dominant eigenvalue in each 
asymptotic sector: 9 

9 This is exactly the same empirical behavior found for the limiting curves B m associated to the chromatic 
zeros of the family 5* m) „ when n^oo [43, Conjeture 7.1]. The graph S TO; „ can be regarded as a square-lattice 
grid with m columns, n rows, free boundary conditions in the longitudinal direction, and special boundary 
conditions in the transverse direction: we introduce two extra vertices such that every vertex on the leftmost 
column of the grid is connected to one of them, and every vertex on the rightmost column is connected to 
the other one. 



30 



Conjecture 7.1 Fix k > 1. T/ien, £/ie limiting curve B k has exactly 2k outward branches 
extending to Q = oo, with asymptotic angles axgQ = 9 k j where 

Okj = p forj = 1,... ,2k. (7.1) 

Moreover, the dominant eigenvalue is (J,k+i,o,* ^ n the asymptotic regions 

a s n^a f S i = l,3,...,2fc- 1 if k is even 

I j = 2, 4, . . . , 2k if k is odd 

while in the other asymptotic sectors the dominant eigenvalue is /ifc+i,i,(i),*- 

We may conjecture that, as k — > oo, the limiting curves B k (without the outward 
branches) converge to some curve B^. This is what Figure 6 suggests. Indeed, outside 
this curve, the outward branches seem to get denser as k increases. Therefore, we conjecture 
that the set of accumulation points of the flow-polynomial zeros for the generalised Petersen 
graphs Gink, k) is dense in the whole complex Q-plane, except in the interior of the curve 
Boo. If this conjecture is true, it implies in particular that flow roots of cubic graphs are 
dense in the complex plane, except possibly in some bounded region near the origin. Note 
that a similar result holds (by duality) for the flow roots of planar graphs [48] . 

If we restrict to the cases with k > 3 (see Figures 3-5), we notice some regularities about 
the isolated limiting points and the dominant eigenvalues in the non-asymptotic region of 
the complex g-plane: 

Conjecture 7.2 Fix k > 4. Then, the dominant eigenvalue in the regions that intersect the 
real Q-axis is 

A*fc+i,i,(i),* for Re Q E (-oo, 2] 
A*k+i,2,(2),* for Re Q e [2,4] 
/'* = < A*fc+i,3,(3),* for Re Q G [4, Q c (k)] (7.3) 
A*fc+i,i,(i),* f° r ReQ e [Q c (k), oo) and even k 
/-ffc+1,0,* f or R e Q [Qc(k), oo) and odd k 

Therefore, Q = 1, Q = 3 and Q = 5 (only when Q c (k) > 5) are isolated limiting points; and 
Q = 2, Q = 4, and Q = Q c (k) are non-isolated limiting points. 

Notice that the above conjecture is also valid for k — 3, except that in this case, there 
is no region where the eigenvalue //fc+i,3,(3),* becomes dominant, as Q c (3) < 4. It is worth 
noting that the dominant eigenvalue comes from the completely symmetric irreducible rep- 
resentation of Se for £ = 1,2, 3. 

Assuming that the dominant eigenvalue on the real Q-axis invariably comes from the 
completely symmetric representation, it is possible to extend the numerical determination 
of Q c (k) to higher values of k. The relevant transfer matrices have dimension (^!) -1 |^fc+il 
and can be diagonalised by a standard iterative scheme that uses their decomposition as a 
product of sparse matrices, as in (3.8). The resulting values of Q c (k) are shown in Table 1. 
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k 


Qc(k) 


1 


3 


2 


3.6180339887 


3 


3.7818423129 


4 


4.5697435537 


5 


4.9029018077 


6 


5.1079785012 


7 


5.2352605291 


o 
O 


o.oz4oyooyUo 


9 


5.3886186958 


10 


5.4364766073 


11 


5.4729804532 



Table 1: Values of Q c {k) for the generalised Petersen graphs G(nk, k) for 1 < k < 11. 



The values of Q c (k) are seen to be well fitted by a second-degree polynomial in 1/k, 
provided that data for k even and k odd are fitted separately. Both fits are consistent with 
each other and give the common limit 

Q c = KmQ c (k) = 5.75 ±0.01. (7.4) 

fc— s-oo 

Formalising slightly these experimental results, we conjecture that the exists a real number 
Q c ~ 5.75 that is an accumulation point of real flow zeros for the graph family G(nk, k), in 
the limit lim^oo lim^^oo, and that no generalised Petersen graph has a real flow zero above 
Q c . The hope that no graph family fares better than the generalised Petersen graph is our 
main motivation — and evidence — for stating Conjecture 1.9. 
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Zeros Petersen k 



Zeros Petersen k = 1 



C\2 



or 



I 



I 




Of 



I 



I 



Zeros Petersen k = 2 



-0.5 

Re(l/Q) 

(b) 

Zeros Petersen k = 2 
l — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — r 




Figure 2: Zeros of the flow polynomial and limiting curves Bt for the generalised Petersen 
graphs G(nk,k) with k — 1 (a,b), and k = 2 (c,d). In (a,c) we show these results in the 
complex Q-plane, while in (b,d), we use the complex 1/Q-plane. For each k, we show the 
flow-polynomial zeros of the generalised Petersen graphs G(12k,k) (□ black), G(16k,k) (o 
red), G(20k,k) (A green), and G(2Ak,k) (0 blue). 



33 



Zeros Petersen k = 3 Zeros Petersen k = 3 




Figure 3: Zeros of the flow polynomial and limiting curves Bt for the generalised Petersen 
graphs P(nk,k) with k = 3 (a,b), and k = 4 (c,d). In (a,c) we show these results in the 
complex Q-plane, while in (b,d), we use the complex l/Q-plane. For each k, we show the 
flow-polynomial zeros of the generalised Petersen graphs P(12k,k) (□ black), P(16k,k) (o 
red), P(20k,k) (A green), and P(24k,k) (0 blue). 
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Zeros Petersen k = 5 Zeros Petersen k = 5 




Figure 4: Zeros of the flow polynomial and limiting curves Bt for the generalised Petersen 
graphs P(nk,k) with k = 5 (a,b), and k = 6 (c,d). In (a,c) we show these results in the 
complex Q-plane, while in (b,d), we use the complex l/Q-plane. For each k, we show the 
flow-polynomial zeros of the generalised Petersen graphs P(12k,k) (□ black), P(16k,k) (o 
red), P(20k,k) (A green), and P(24k,k) (0 blue). 
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Figure 5: Zeros of the flow polynomial and limiting curves Bj for the generalised Petersen 
graphs P(7n, 7). In (a) we show these results in the complex Q-plane, while in (b), we 
use the complex 1 / Q-plane. We show the flow-polynomial zeros of the generalised Petersen 
graphs P(21,7) (□ black), P(35, 7) (o red), P(63, 7) (A green), and P(119,7) (0 blue). 



8 Discussion 

In this concluding section we would like to place our results on the flow polynomial in 
a larger perspective, and to propose a number of conjectures that should stimulate further 
research. We have seen that the Q-state Potts model with temperature parameter v com- 
prises both the flow polynomial (v = —Q) and the chromatic polynomial (v = —1) as special 
cases. It is therefore useful to compare our results on the flow polynomial for the particular 
class of non-planar graphs G(n, k) with what is known for the chromatic polynomial for 
certain classes of planar graphs. More generally, we should like to outline and compare the 
results and expectations for real accumulation points of chromatic zeros for planar graphs, 
and of flow zeros for non-planar graphs, when the graphs are taken to consist of n identical 
layers of width L vertices and we impose periodic boundary conditions in the n-direction. 
Going beyond this, we discuss also the behaviour of real accumulation points of zeros for the 
Potts-model partition function Za(Q,v) along more general curves in the (Q,v) plane. 

Since many of the recent results on chromatic zeros have been obtained at the borderline 
between mathematics and theoretical physics, we shall use — for lack of a better term — 
the denomination "result" to characterise a fact which has been derived using exact (but 
non-rigorous in the mathematical sense) methods of theoretical physics. These methods 
include in particular two-dimensional conformal field theory and integrability techniques, 
and renormalisation group arguments. 

Our main objective is to spell out a number of striking analogies between, on the one hand, 
the intersection of the limiting curves of Figure 6 with the real Q-axis, and Conjecture 7.2 
on the dominance pattern of eigenvalues on the real Q-axis, and, on the other hand, the 
corresponding results [26] for the chromatic zeros for cyclic strips of regular planar lattices. 
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Limiting curves Petersen Limiting curves Petersen 




Figure 6: Limiting curves for the generalised Petersen graphs P(kn, k) in the limit n — > oo. 
In (a) we show these results in the complex Q-plane, while in (b), we use the complex 1/Q- 
plane. We show the curves corresponding to k — 1 (black), k = 2 (red), k = 3 (green), k = 4 
(blue), k = 5 (magenta), k = 6 (orange), and k = 7 (gray). 



In the planar chromatic case the general situation is expected to be the following: 10 

Conjecture 8.1 Consider a family G n ^ of planar connected graphs, consisting ofn identical 
layers of width L vertices, with periodic boundary conditions in the n-direction. Let Al (resp. 
Bl) denote the set of limiting points of chromatic zeros for G n ^, in the limit n — > oo, that 
satisfies condition (a) (resp. condition (b)) of Theorem 5.1. Further let B p be the p-th Beraha 
number: 

B p = 4 cos 2 Q , peN, (8.1) 

and let 

B cvcn = {Q = Bp : pe2N} , B odd = {Q = B p : p G 2N + 1} 
denote, respectively, the set of even and odd Beraha numbers B p . Then: 

1. There exists, for each L, a number Q C (L) e (0,4) so that 

(a) A L f]R = B odd n [0, Q C (L)) 

(b) B L n R = £? cvcn n [0, Q C (L)) U {Q C (L)} 

2. For any p & N , the dominant eigenvalue of the transfer matrix for 

Q e (B 2p ,B 2(p+1) )n [o,Q c (L)) 

comes from the sector with p marked clusters. 
10 In this section, B v denotes the Beraha numbers (8.1), not the Bell numbers B n (3.10) defined earlier. 
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Compelling evidence for Conjecture 8.1 comes from the 

Proposition 8.2 (Jacobsen-Salas [26]) Conjecture 8.1 is correct for cyclic strips of the 
square and triangular lattices with L < 8. The corresponding values of Q C (L) can be found 
in [26]. 

From the point of view of statistical physics, the natural step suggested by these finite- L 
results is to go to the so-called infinite-volume limit, lim^oo. We have then the following 
non-rigorous "results" obtained by the methods of theoretical physics: 

Result 8.3 [26,44,45] For each family of planar graphs G Ht L, defined as above, the limit 
Q c = lim^oo Q C (L) exists and satisfies < Q c < 4. 

Result 8.4 [4—6,28,29] For cyclic strips of the square lattice, Q c = 3. For cyclic strips 
of the triangular lattice, Q c = 2 + \/3. 

The deep meaning of these results is that the special role of the Beraha numbers in the 
production of real chromatic zeros is brought out whenever the chromatic line (v = —1) 
intersects the so-called Berker-Kadanoff (BK) phase in the (Q, v) plane of the Potts model 
[44,45]. For each graph family, and for each Q e [0, 4), it can be argued [30] that the extent of 
the BK phase is a finite non-empty interval v_(Q) < v < v + (Q) with the properties f+(0) = 
and limQ^4f + ((5) = liniQ-^ v -(Q)- In particular, other choices of the temperature curve 
v(Q) are possible and might in some cases lead to higher values of Q c than those in Result 8.4. 
The bound on Q c given in Result 8.3 is even sharp: 

Result 8.5 There exists a family of planar graphs G n ^, defined as above, and a curve 
v = f(Q) in the (Q,v) plane, so that the limiting set of zeros for the Potts-model partition 
function Z G {Q,v) satisfies Conjecture 8.1 and Result 8.3 with Q c = 4. 

Indeed, this is the case [30,44,45] for cyclic strips of the square lattice with the choice [2] 
v 2 = Q for < Q < 4 and —2 < v < 0. Another example is provided by cyclic strips of the 
triangular lattice with the choice [3] v 3 + 3v 2 = Q for < Q < 4 and — 2 < v < 0. 

Let us finally mention that Conjecture 8.1 enjoys strong support from the special repre- 
sentation theory of the quantum group at roots of unity [44,45]. 

Turning now to the flow zeros of non-planar graphs, Figure 6 and Conjecture 7.2 makes 
us confident that the statements of Conjecture 8.1 hold true also in the flow case, provided 
that one replaces even/odd Beraha numbers by even/odd integers. On a fundamental level, 
the replacement of Beraha numbers by integers stands out most clearly by comparing the 
eigenvalue amplitudes in the two cases [see, e.g., [26, Eqs. (2.28)-(2.29)] and Eq. (3.14)]. To 
be precise: 

Conjecture 8.6 Consider a family G n ^ of non-planar connected graphs, consisting of n 
identical layers of width L vertices, with periodic boundary conditions in the n-direction. Let 
Al (resp. El) denote the set of limiting points of flow zeros for G n ^, in the limit n — > oo, 
that satisfies condition (a) (resp. condition (b)) of Theorem 5.1. Then: 

1. There exists, for each L, a number Q C (L) > so that 
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(a) Al fi IR = (2N — 1) n [0, Q C (L)) 

(b) B L n K = (2N) n [0, Q C (L)) U {Q C (L)} 

2. For any p G N, the dominant eigenvalue of the transfer matrix for 

q e (2(p-i),2 P )n[o,g c (L)) 

comes from the fully symmetric irreducible representation (p) of the sector with p 
marked clusters. 

This conjecture has been validated in this paper for G n ^ = G(nk,k), the generalised 
Petersen graphs, with L = k + 1 and k < 7. To parallel Results 8.3-8.4, we have presented 
numerical evidence that the limit Q c exists in this case, with the value (7.4). 

Improving on Conjecture 8.6 we can characterise more fully the structure of the eigen- 
values and their crossings. 

Conjecture 8.7 Consider a family G n ^ of non-planar connected graphs, as above. Let Qi 
denote the smallest Q > so that the transfer matrix T — corresponding to the fully 
symmetric representation (£) of £ marked blocks — has two or more dominant eigenvalues, or 
let Qe = oo if no such Q exists. For Q e (0,Qe), let \^ denote the corresponding (unique) 
dominant eigenvalue. Then: 

1. Q < Qi < Q 2 < ■ ■ ■ ■ 

2. \e is real. 

3. The equation |A^| = |A^| has at most one solution for < Q < mm(Qe, Qe). 
4- If such a solution exists, it is Q — £ + £' — 1. 

This conjecture has been validated in this paper for G n i = Gink, k), with L — k + 1 and 
k < 7. 

There is solid evidence [6, 28, 41] that the planar graphs with the largest possible real 
chromatic roots are large pieces of the triangular lattice. However, when cyclic boundary 
conditions are imposed, even such lattices cannot saturate the bound in Result 8.3 (i.e., have 
Q c — 4), as we have seen in Result 8.4. Replacing now chromatic roots of planar graphs with 
flow roots of non-planar graphs, a comparison between Conjectures 8.1 and 8.6 reveals that 
the special role of Q = 4 in the former case is replaced by Q = oo in the latter. But although 
we have studied the generalised Petersen graphs because they could be expected — on the 
basis of sporadic evidence for small graphs [22] — to have the largest possible real flow roots, 
we should not forget that the flow polynomial is merely another special case of the Potts 
model. Therefore, by analogy, there seems no reason to believe that real flow roots can 
be arbitrarily large. This motivates, once more, Conjecture 1.9 that we have made in the 
introduction. 

To recover the complete analogy with Result 8.5 we should therefore like to state one last 

Conjecture 8.8 There exists a family of non-planar graphs G Hj l, defined as above, and a 
curve v = f(Q) in the (Q,v) plane, so that the limiting set of zeros for the Potts-model 
partition function Z G (Q,v) satisfies Conjecture 8.1, but with lim^oo Q C (L) = oo. 
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To conclude, we recall that for the Potts model on planar graphs G n: L, when Q = 
(2 cos(7r/p)) 2 is a Beraha number the representation theory of the underlying quantum alge- 
bra [38] leads to massive degeneracies in the spectrum of the transfer matrix. In particular, 
depending on p, complete sectors of eigenvalues are contained in other sectors in an inclusion- 
exclusion fashion, completely independent of the value of v. 

It seems quite plausible that something similar should happen in the non-planar case, 
when Q is an integer. We have indeed observed such phenomena for the transfer matrix 
(3.8) of the Potts model on the graphs G(n,k). Let us mention but the simplest example 
of our findings, which happens for Q = and any v. The set of eigenvalues in the £ = 
sector is then fully contained in the set of £ = 1 eigenvalues. Moreover, let Wg denote the 
set of eigenvalues in the fully antisymmetric representation (1, 1, . . . , 1) of the sector £ + 1. 
We have then observed, for any I > 1, that Wg+i \ Wg +2 C W t . 

A Proof of the lemmas related to the transfer matrices 

In this appendix we will provide the proofs of the extensions of Lemma 4.1, Lemma 4.2, 
and Corollary 4.3 for all k > 1. We will also prove Lemma 4.6. 

It is useful to rewrite the transfer matrix T L (3.8) as the product of two operators 
T L = H V, given by 



H = n H o* V ° = H 01 VoHo 2 V H 03 ---V H ofe Vo (A.la) 

i=k 
k 

V = J] V * = ViV 2 ...V fc _!V fc (A.lb) 



i=i 



where L — k + 1, the operators V, and H 0i are defined in (3.7) with v = —Q, and we have 
used the property [V,, V^] = [Dj, Dj] = for all 

A.l Extension of Lemma 4.1 

Lemma A.l Fix k > 1, and I — 0, . . . , k + 1. Then, the diagonal block Tk+1,1 contains the 
trivial eigenvalue Hk,k+i = ( — l) fe with multiplicity Nk,\{£)£\, where: 

'0 if£ = 

k 



^■j^^-i-sOfc-i}^ (A,2) 

and S n is given in (3.33). 



PROOF. Let us consider first the transfer matrix Tl — PzTl with Tl given by (A.l), and the 
projector defined in (3.17). Let us fix £ > 1 and a representation A G Si of dimensionality 
dim A. The set can be split into two disjoint sets: 
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• C[ is the set of all partitions of {0,1, ... ,k} with no un-marked singletons, £ marked 
clusters, and vertex is a marked singleton. 

• C2 is the set of all partitions of {0,1,..., k} with no un-marked singletons, £ marked 
clusters, and vertex is not a singleton. In this case, the vertex can be either marked 
or unmarked. 

It is clear that \A^\ = \C[\ + . Notice that we do not need to consider the case of vertex 
being an un-marked singleton, as all partitions of this type do not belong to 

In order to find out the structure of the transfer matrix T^, we do not need to specify 
the representation A; therefore, we will omit this index (as well as the index k) to make the 
notation clearer. 

Let us first consider a partition V\ G C[ represented by the vector e-p L in the space of 
partitions. Then, 

Me Vl = {-Q) k e Vl +J2 A -r e Vi ( A - 3 ) 

Vecf 

where the sum is over the set of all partitions C[ of {0, 1, . . . , k} with £ marked blocks, with 
vertex a marked singleton, and such that there is at least one un-marked singleton. This 
is because V only contains the identity operator and detach operators Dj with 1 < i < k. 
If we consider a partition P 2 G C|j then 

V e V2 = {-Q) k e V2 +^ A v e v +^ B v e v , (A.4) 

reef va& 2 

where now is the set of all partitions of {0, 1, . . . , k} with £ marked blocks, with vertex 
not a singleton, and such that there is at least one un-marked singleton. We omit the 
possible partition with an un-marked singleton, as it belongs to the kernel of the first 
operator V in H. 

We now notice that the operator H jV can be written as: 

H 0l V = -Q/ + Do + Q 2 J 0l -QJoiD (A.5) 

where Joi and D are the join and detach operators, respectively. It is not hard to see that 
for all partitions V G A^\ 

Y\ 0i \l e v = -Qe P + ^ A v , e v > , (A. 6) 

V'ecf 2 

where if V G C^, V might coincide with V . For partitions with at least one un-marked 
singleton V G Cf U C^, we have 

HoA/oe-p = -Qe v + ^ A v > e v * + ^ B v < e v > . (A. 7) 

vec l 2 v'ec f 2 

Notice that when we apply the full operator H, the partitions belonging to C[ U are 
eliminated by the application of the operator P L (3.17). 
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We now put together the above observations, and conclude that if V e C[, then 

J L e P = {-Q) 2k e v + A P , e v , . (A.8) 

7"ec| 

On the other hand, if V G C l 2 , then 

T L e P = ^ Ap/ ep/ . (A.9) 

This means that, if we order the partition states appropriately, the transfer matrix T l has 
a block triangular form. This result holds for all representations A e 5V 

The block corresponding to the partitions in C[ is diagonal with all eigenvalues equal to 
Q 2k . In terms of the matrix T L (3.20), this common eigenvalue takes the value fik,k+i — 
(—l) k Q~ 2k Q 2k = (— l) k . This proves the existence of the trivial eigenvalue Hk,k+i = {— 
off L . 

The dimension of this diagonal block should be equal to the number of partitions of the 
set {1, 2, ... , k} with no un-marked singletons, and £ — 1 marked blocks. (Indeed, for £ = 0, 
this dimension is identically zero.) But now, as we are considering linear combinations 
of the partitions with the right symmetries under Si, these blocks should be considered 
indistinguishable. Therefore, 

\Ci\ = ^^ir^dimA) 2 , (A.10) 

as there are dim A Young tableaux of dimensionality dim A. For a fixed value of £ > 1, the 
total number of trivial eigenvalues in Tl coming from all irreducible representations A e Si 
and divided by £\ is 

w = I. IF^Ty.^ g dlm ' 2 A = jrhy:*"'' • (A ' n) 

as ^xeSe dim 2 A = i\. Finally, using (3.34) and after some algebra, it is not difficult to find 
the expression (A. 2). I 



A. 2 Extension of Lemma 4.2 

Lemma A. 2 Fix k > 1, and £ = k. Then, there are k non-trivial eigenvalues Hk,k,s ^ n the 
diagonal block T fc fc ( once we have removed the trivial eigenvalue Hk,k+i = {-If). Each of 
these eigenvalues Hk,k,s has multiplicity k\. 

Proof. Let us first consider the following partition (010'1')(22 / ) . . . [kk') with £ = k links, 
where the primed numbers denote the sites on the top row, and the un-primed ones, those 
on the bottom row. This partition does not lead to the trivial eigenvalue fik,k+i — ( — ^) k , as 
it belongs to the class C k defined in the proof of Lemma A.l above. Then the action of the 
transfer matrix T^+i [cf., (A.l)] on this partition generates exactly k partitions that can be 
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constructed as follows: starting from the basic partition (Oil') (22') . . . (kk') with £ = k links, 
attach the site 0' to any of the k possible blocks. It is important to note that the ordering of 
the k links is preserved in the whole process. (Loosely speaking, there is "no room" to switch 
two links if we have k + 1 sites and k links.) This also means that the action of the transfer 
matrix T k +i on any of these partitions is independent of the actual ordering of the links. 
Hence, for each of the k\ possible orderings of the k links, the transfer matrix will be the 
same (modulo some reordering of the partitions), and will have dimensionality k. Therefore, 
all irreducible representations of the symmetric group S k will give the same k non-trivial 
eigenvalues, and the multiplicity of each of these eigenvalues in T fc+l fe is Xaes dim 2 \ — k\. 



A. 3 Extension of Corollary 4.3 

Corollary A. 3 Fix k > 1. Then for all pairs (I, A) with < £ < k — 1 and A e S^, the 
diagonal block T^+i^a has N kj0 (£) non-trivial eigenvalues [i k ,e,\,s (once we have removed the 
trivial eigenvalue fik+i,k = (— l) fe ; which appears for every 1 < £ < k + 1). The value of 
Nk,o(£) is given by: 



p=0 ^ ' 



N k , (£) 



if£ = 



P =i 



k\ \v 

P 



[S k - P + S k+1 _ p ] if\<£<k-\ 
if£ = k 

lf£=k+l 



(A.12) 



where S n is given by (3.33), and Yg is the number of Young tableaux with £ cells. 



Proof. Let us first fix < £ < k, and the irreducible representation A of the symmetric 
group Se (of dimensionality dim A). The corresponding diagonal block is T k +i,e,\, and its 
dimensionality is 

dim W = ^ M ' ■ 



£\ 



dim A . 



(A.13) 



The number of trivial eigenvalues /i ktk +i = (— l) k in T^+i^ is N kt i(£) dimA by Lemma A.l. 
Therefore, the number of non-trivial eigenvalues in this block is given by 



N, 



£\ 



- N Kl {£) 



dim A , 



(A.14) 



and the number of non-trivial eigenvalues for all different values of A G St is given by 

~N k {£) 



N k , (£) 



£\ 



N k:1 (£) 



dim A . 

\es e 



(A.15) 
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This latter sum is equal to the number of Young tableaux with £ cells Yg (4.6). 

We can provide a closed form for N k:0 (£) for < £ < k by combining Lemma A.l, the 
definition N k {£) = \A k e) +1 \, and (3.34)/(A.15). After some algebra we find Eq. (A.12) for 
£ < k. 

The case £ = k is derived directly from Lemma A. 2. There are k non-trivial eigenvalues 
for every irreducible representation A of S k . Therefore, N kj0 (k) = k. 

Finally, for £ — k + 1, the number of non-trivial eigenvalues is zero, as all eigenvalues are 
trivial in this sector. Note that N k (k + 1) = (k + 1)! [cf. (3.34)], and that N k ,i(k + 1) = 1 
[cf. (4.4)], so that N kfl (k + 1) = by Eq. (A. 15). I 

A. 4 Proof of Lemma 4.6 

All the matrix elements of T k +i t e,\ are polynomials in Q. Therefore, tr T^ +1 e A will be a 
polynomial in Q. However, not all the matrix elements of T^+i^a are polynomials in Q; so 
we cannot conclude that tr T^ +1 e A is a polynomial in Q. 

We need a different argument. Let us consider any partition V belonging to A^ k+1 . Then, 

its contribution to trT k+1 ^ x will be the sum of the contributions of all the diagrams that 
start with the partition V and end with the same partition after n steps. By inspection of 
the form of the transfer matrix T k +\ (3.8) and its components Vj and Hoi (3.7), it is clear 
that there is a contribution proportional to {—Q) 2kn due to the application of all the identity 
operators in T^+i (for each layer, there are 3k of them, and only 2k take the factor —Q). 
Indeed, this is the minimum number of — Q factors one can possibly obtain for a diagram 
of this type. Therefore, the minimum power of Q that appears in trT^ +i e x is Q 2kn , which 

is exactly the prefactor in the definition of T fc+i (3.19)/(3.20). In conclusion, tiJ k+iex is 
indeed a polynomial in Q. 

Once the polynomial character of tr T^ +1 £ A is established, we have to take care of the 
degree of this polynomial. The lemma is proved if we are able to show that all the entries 
in Tk+i,e,x have powers of Q of degree at most 3k + min(l — £,0). This would imply, that 
trT^ +1£ A is a polynomial in Q of degree at most n[3k + min(l — £,0)}, and that trT£ +1£A is 
a polynomial in Q of degree at most n[k + min(l — £, 0)], as claimed. 

Let us first consider the case £ = 0. Let us start with an arbitrary partition of the top 
row V G ^4.^+!- When we apply T^ +1 to that partition, the maximum number of — Q factors 
that can appear is (— Q) 3k : we apply the part — Ql for each vertical operator Vj, and the 
join operator — QSoi for each horizontal operator Hoj. We obtain for all V G ^.j^ the same 
final partition: {{0, 1, . . . , k}}. Therefore, tr T k+1 A is a polynomial in Q of degree at most 
3kn. 

The case £ = 1 is similar. If we start with an arbitrary partition of the top row V G A^ 1 ^ 
with marked, and we apply the same operators as above, we end up with the partition 
{{0,1, ... ,k}} (where the overline over a site denotes that this site is connected to the 
bottom row). The coefficient corresponding to this partition is a polynomial in Q of degree 
at most 3k. Therefore, trT k+l l x is a polynomial in Q of degree at most 3kn. 

The case 2<£<A; + lis similar. Let us start with the following simple partition in A k + ± : 
{{0,1, . . . , k}, {!}, ...,{£— 1}. The argument is similar to the previous cases, except that 
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we cannot join the site to any of the other £ — 1 blocks; since if we did, then the number 
of links would be smaller than £. Therefore, we get a diagonal entry which is a polynomial 
in Q of degree at most 3k — (£ — 1). Therefore, trT k+iex is a polynomial in Q of degree at 
most n[3k- (£ - 1)]. I 

Remark. Note that the above Lemma for k = £ + 1 implies that the trace tr T^ +l k+l x is a 

polynomial in Q of degree at most n[3k — k] — Ink. Thus, tr T£ +1 A is a polynomial in Q 
of degree at most 2nk — 2nk = 0, in agreement with Lemma A.l. 

B The polynomial 3>g(ii9,7) 

_J Oj 

240453758183717079931230416441214627161100181583221695778758847017660 

1 3778010581676303383947166862404894626185168386864045862610150115812052 

2 30335569899732630785396756910315613411372095140384361997526113182374163 

3 165825101972051263346220423069919998872239305113596582973635632055394977 

4 693579540994844285783536577772823633624035569377766733899721814879025261 

5 2365298420910039361031441041999881902182643546761840393901756834507719194 

6 6843881667041711337123509683076504016112680686224789087846006938716709936 

7 17263457862032800683223779137831846913199242114984324936441639648268053289 

8 38713433713214705146001095156944656719927970414458008688062316276391488388 

9 78323139898934533680532026518046626525971702819240021820160355852376101900 

10 144596924706446634397013884950643553134413962554815223156161387400333508225 

11 245805490880525223783583677086377232647690995050227537846281841251646914824 

12 387578029394911704914155641300482722338098395075782035706404804765273096365 

13 570251442350241862197997781738305884386349966454124613787648726411773264858 

14 786830079084979854704322218255681568053263515605219068314256695784834361687 

15 1022415665287852311682199189937274383849816815595817077282491287814210089327 

16 1255617405561053292049115611630960519206807161888631329237623843787650902598 

17 1461844034312439449575494091690068905133069052914279852509149780627087441657 

18 1617740821249798765950531657277432254427341252461560806920840933023592630423 

19 1705621288663026987807981694789591864355698512752538498113371227662898599836 

20 1716721077841948079359918580031477807370750494656542267873960267539031216693 

21 1652471749735377739204603329938965183162646296909313884512941374410264086521 

22 1523600861689858982867164652432211354402872413669859092087649541160917763276 

23 1347478648579683679991904590350043152972817068987716895096810710111624258414 

24 1144543577841773436389066666439373326178625474359360165078175860151437888418 

25 934744903425435172192816740334988586706047594412696933301815485640045053088 

26 734761463017109267283189805699271455927373168195072916098687225728186938239 

27 556406139434285183296690956661070857540502747488343675526234531927677714022 



Table 2: Coefficients of the flow polynomial $G(ii9,7) (Q) = (Q ~ 1)(Q ~ 2 )(Q - 

117 

3)£(-l) i+ W 

i=0 
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J Oj 

28 406248919977322871211475078592598791327840251152453917692944363153498710100 

29 286204685555683178005609409552066982429730669340777755008844846916076522994 

30 194690835588052864277724799113473665492499996408415866330514336755095511519 

31 127959409173717651616971568749651845945079328780268936837379415435944126802 

32 81303058361075417592304554906219204545070829443708169095378279034181616797 

33 49966395799032436266189002396534646533235474102922584729978792549768709570 



34 29716297381875173070101426039892101696811371497930415707024934873184994971 

35 17109926011423464771467167427507814898291755276949983853151230858292134005 

36 9541424411645170221820026129682497673840857753028326218052776323947917867 

37 5155250259203698010398273841766642027281260863127250587163264787844562005 

38 2699625968148386503782707617994694138490854643501950268752229332983496484 

39 1370590172276847791788823798202187745155060493411059843973697377813429115 

40 674811604481186691941181395279228274700479647745836519066846509110686871 

41 322283263195622661124593258235203925484545860053525952388302407768566612 

42 149338540704756141327132733140880585247255726988170348306522536601036428 

43 67154254327764896250779574488797365684816293909986077788680077649219059 

44 29310387712861712336349988437753297027267572645934195793910184667409525 

45 12418975519744799441679142496250486367091770751393334192303593690937216 

46 5108897246470415157235242805025854028977310258147102091244333174463070 

47 2040790191809533248688898356534425242558953177861065156649563146370679 

48 791669367678994305215119572575588439395174310193863932545018547493949 

49 298263611514781343272224652335166519468558609433860386772388696144855 

50 109143202315633324987330964029233262735802381917881292972386146290583 

51 38793120529207946517836936596295050010421179961661603601012761033976 

52 13393333822523527127943242411311403493467668611668190366289160119798 

53 4491649304261905729281333948467503543032844169847443690259485449787 

54 1463208911995358578670442816667443146618570037914577188051167519252 

55 463002576567921398643877003491825254798116352410207026582354227639 

56 142306049793794840801481088729855032253716088115496415256458286993 

57 42482066138212698853729014913405235508231962094892775535623301954 

58 123 1693782063003445771766 13 1 18 13898450293797908401593794293724758 

59 3468010035454592592294532401090077198723295654493559716526595085 

60 948193746090263874793189178585310649056926723273402296835255341 

61 251712395999861219430022121008687006528308619159639108633408792 

62 64870522436154805702181162937007326080100049608985554185883024 

63 16227914273399692643538894326321813038077624088578343763024693 

64 3939856175099340812778271579647788724054817562072496203271622 

65 928159554704604760538186427442530041159446537967362956747430 

66 212130925228863323822024721431112601866613699451131860235710 

67 47025221622368034908085340305586428332317539464533841648123 

68 10108841512034189983877166351043431883332647922883010856327 

69 2106711947665769238334508701464763203793865636828903442294 

70 425522761043753110078251752857561183613194192306420512412 

71 83277105914548332097899604836384635641176843712157856926 

72 15786114232928385188569755800951974220820338221822797673 

73 2897497349085917262770771233534175446735608061805601330 



Table 2: (Continued.) 



_j_ Oj, 

74 514765848404059046836837761417162573532059958487097353 

75 88483984104925516580422896030293322509203586666446608 

76 14709751385762293879136933265239581784434168823726105 

77 2363940443321701747333111550458096052434867769928216 

78 367070636910713149669510355126240632283426228583990 

79 55045491030169621973045797552387375721405805002210 

80 7967363593383576742789482483060195429602666353792 

81 1112432655029526956212696750602542449173138786358 

82 149736221917777447936668185649515326795606657139 

83 19417037687233804941366229781880098811952696949 

84 2423988165427613282917222688678473423362471885 

85 291095664115059080227254369870423630080267269 

86 33600127040413356951575245326064861734434090 

87 3724448845910258224439985334986604271814630 

88 396083067053161605214005007534380917557302 

89 40370994515077502586895936288134103918883 

90 3939422436008602219258246065351204112198 

91 367585380265857178359590490307505526698 

92 32755670407906919539877604731444999292 

93 2783619278733903044429223967432069173 

94 22525141385362067800441 1 18 1 186579261 

95 17327620037637857872926704478452708 

96 1264839938089389154908070831758552 

97 87436328017906227077510426235281 

98 5711550730808768130994711119064 

99 351695047292333601570905687819 

100 20358915651379087237820736779 

101 1104604810216428455300363966 

102 55981869638360027493471673 

103 2640008664516116400095478 

104 115338644735668374521021 

105 4644751096461226418781 

106 171402937526338140238 

107 5756292666737489361 

108 174485239351077202 

109 4726362454421551 

110 113000004047458 

111 2347463061581 

112 41510490693 

113 607499109 

114 7064107 

115 61191 

116 351 

117 1 
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